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(57) Abstract 

A method and apparatus for the exploration of hydrocarbons comprising the steps of: obtaining a set of seismic signal traces distributed 
over a predetermined three-dimensional volume of the earth; dividing the three-dimensional volume into a plurality of analysis cells having 
portions of at least two seismic traces located therein; computing outer products of the seismic traces within each cell; forming the covariance 
matrix for each cell from these outer products; computing the dominant eigenvalue and the sum of the eigenvalues of the covariance matrix 
of each cell and computing a seismic attribute from the ratio of the dominant eigenvalue to the sum of the eigenvalues of the covariance 
matrix of each cell; and forming map of the seismic attributes of selected groups of cells. 
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Method and Apparatus for 
SEISMIC Siqmal Processing and Exploration 

Technical Field 

5 This invention relates to the general subject of seismic exploration and, in 

particular, to apparatus and methods for. the exploration and production of oil and 
gas by identifying structural and stratigraphic features in three dimensions. 

Background Of The Invention 

10 in seismic exploration, seismic data is acquired along lines that consist of 

geophone arrays onshore or hydrophone streamer traverses offshore. Geophones 
and hydrophones act as sensors to receive energy that is transmitted into the 
ground and reflected back to the surface from subsurface rock interfaces. Energy 
is usually provided onshore by Vibroseis® vehicles which transmit pulses by 

15 shaking the ground at predetermined intervals and frequencies on the surface. 
Offshore, airgun sources are usually used. Subtle changes in the energy returned 
to surface often reflect variations in the stratigraphic, structural and fluid contents of 
the reservoirs. 

In performing three-dimensional (3-0) seismic exploration, the principle is 
20 similar however, lines and arrays are more closely spaced to provide more detailed 
sub-surface coverage. With this high density coverage, extremely large volumes of 
digital data need to be recorded, stored and processed before final interpretation 
can be made. Processing requires extensive computer resources and complex 
software to enhance the signal received from the subsurface and to mute 
25 accompanying noise which masks the signal. 

Once the data is processed, geophysical staff compile and Interpret the 3-D 
seismic information in the form of a 3-D data cube (See FIG. 1) which effectively 
represents a display of subsurface features. Using this data cube, information can 
be displayed in various forms. Horizontal time slice maps can be made at selected 
30 depths (See FIG 2). Using a computer workstation an interpreter can also slice 
through the field to investigate reservoir issues at different seismic horizons. 
Vertical slices or sections can also be made in any direction using seismic or weil 
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data. Time maps can be converted to depth to provide a structural interpretation at 
a specific level 

Seismic data has been traditionally acquired and processed for the purpose 
of imaging seismic reflections. However, changes in stratigraphy are often difficult 
5 to detect on traditional seismic displays due to the limited amount of information 
that stratigraphic features present in a cross-section view. Although such views 
provide an opportunity to see a much larger portion of these features, it is difficult to 
identify fault surfaces within a 3-D volume where no fault reflections have been 
recorded. 

10 Coherence and semblence (a measure of multichannel coherence) are two 

measures of seismic trace similarity or dissimilarity. As two seismic traces increase 
in coherence, the more they are alike. Assigning a coherence measure on a scale 
from zero to one. "0" indicates the greatest lack of similarity, while a value of "1" 
indicates total or complete similarity (i.e.. two identical traces). Coherence for more 

1 5 than two traces may be defined in a similar way. 

One method for computing coherence was disclosed in a U.S. Patent 
Application by Bahorich and Farmer (assigned to Amoco Corporation) having a 
serial number of 067353,934 and a filing date of December 12, 1994. A method for 
computing semblance was disclosed in a U.S. Patent Application by Marfurt et al. 

20 (assigned to Amoco Corporation) having an application number of 60/005.032 and 
a filing date of October 6. 1995. The Marfurt et al. invention included a brute force 
search over candidate dips and azimuths. 

As good as both methods have proved to be. they have some limitations. 
Improved resolution and computational speed are always desirable. 

25 

Slimr" af Y Qf The Invention 
In accordance with the present invention, a muWtrace eigen-decomposition 
process is disclosed that is more robust and that has higher resolution than 
previously known methods. In one embodiment of the invention, a method is 
30 disclosed for the exploration of gas and oil. The method comprises the steps of: 
accessing a data set of seismic signal traces distributed over a predetermined 
three-dimensional volume of the earth; in a running window determining the outer- 
product of at least two data vectors formed from at toast two seismic traces; 
forming a covariance matrix by adding the outer-products; computing a seismic 
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attribute that is a function" of at least the dominant eigenvalue of the covariance 
matrix; and forming a map from the computed seismic attributes over at least a part 
of the pre-determined three dimensional volume of the earth. 

In another embodiment of the invention, the process of the invention is 
5 encoded on a computer readable media (e.g., magnetic disk, magnetic tape, 
CD-ROM, etc.) for directing the operation of a computer to compute the seismic 
attributes. In other embodiments of the invention, a map is prepared from the above- 
described process and the map is used to locate oil and gas deposits. 

This technique is particularly well suited for interpreting fault planes within a 
10 3-D seismic volume and for detecting subtle stratigraphic features in 3-D. This is 
because seismic traces cut by a fault line generally have a different seismic 
character than seismic traces on either side of the fault Measuring seismic trace 
similarity, (i.e., coherence or 3-D continuity) along a time slice reveals lineaments of 
low coherence along these fault lines. Such coherency values can reveal critical 
15 subsurface details that are not readily apparent on traditional seismic sections. 
Also by calculating coherence along a series of time slices, these fault lineaments 
identify fault planes or surfaces. 

Numerous other advantages and features of the present invention will 
become readily apparent from the following detailed description of the invention, the 
20 embodiments described therein, from the claims, and from the accompanying 
drawings. 



Brief Description of the Drawings 

FIG. 1 is a pictorial representation of the information obtained from 
25 processing 3-D seismic data; 

FIG. 2 is a pictorial representation of a horizontal time slice of 3-D seismic 
data processed in accordance with the prior art 

FIG. 3 is a drawing depicting two adjacent seismic traces; 

FIG's. 4 through 8 are schematic diagrams depicting the coherency of a 
30 pair of seismic traces in accordance with the present invention; 

FIG. 9 is a pictorial representation of a running window analysis cube; 

FIG's. 10A, 10B and 10C are schematic diagrams of groups of seismic 

traces; 
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RG's. 11A t 11B and 11C are schematic diagrams of two dimensional 
analysis windows; 

FIG. 12 is an elementary process flow diagram; and 
FIG's. 13A, 13B and 13C are pictorial representations of the same 
5 horizontal time slice in accordance with the inventions of Bahorich et al., Marfurt et 
ai. and the present invention. 



PfttftiM neserintion 
While this invention is susceptible of embodiment in many different forms, 
10 there is shown in the drawings, and will herein be described in detail, several 
specific embodiments of the invention. It should be understood, however, that the 
present disclosure is to be considered an exemplification of the principles of the 
invention and is not intended to limit the invention to the specific embodiments or 
algorithm so described. 
15 Before describing the process of the invention, the underlying principle of 

operation will be described. 

Consider two traces t, and tj over a specified time window or depth 
window of N samples for which coherence is to be evaluated. A representative 
diagram of the traces and the relevant analysis window is shown in FIG. 3. The 

20 first trace t, consists of the time series (t 11t t 12 t tN ) and the second trace ta 

consists of the time series (t 21f t» W In these two time series, the first index 

refers to the trace number (i.e., trace 1 or trace 2), while the second index refers to 
the sample number. 

By plotting one trace against the other in the familiar two-dimensional 
25 Cartesian coordinate system a better understanding of the meaning of coherence in 
the context of the present invention can be obtained. Plotting equivalent time 

samples from the two traces fl.e.. the point pairs (t 11( t 21 ), (t„, t^) (t 1M , 

produces a crossptot of the two time series. Letting the x-axis represent the first 
trace t, and the y-axis the second trace results in the diagram* depicted in FIG. 
30 4. It is the pattern formed by these points which manifests the coherence of the two 
traces. 

The general shape of these two correlated traces is a set of points 
represented by an ellipse. This ellipse is a generalization since it does not 
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represent each individual point but rather describes the "overall" nature of all the 
points. The major and minor axes of this ellipse will be oriented in a direction which 
is determined by the geometry of the paired points. The lengths of the two axes are 
also determined by this geometry. A typical representation of these points and the 
5 corresponding ellipse is shown in FIG. 5. 

The directions and magnitudes of the major and minor axes of the ellipse 
may be represented by two scaled vectors with the longer vector oriented along the 
major axis and the shorter vector along the minor axis. The magnitudes of these 
two vectors correspond to the two eigenvalues of the data covariance matrix and 
10 the normalized vectors correspond to the eigenvectors. The eigenvectors, scaled 
by their respective eigenvalues, denote the magnitudes and directions of the major 
and minor axes. The "principal component" corresponds to eigenvector which is 
associated with the dominant eigenvalue. 

The figures that follow (FIG's. 6 through 8) are intended to give an intuitive 
15 understanding of the mechanics behind the previous discussion. In these figures, 
traces are constructed using simple amplitude and phase variations, and the effect 
of these variations is observed on the associated eigenvalues and eigenvectors. 

FIG. 6 demonstrates how two identical traces degenerate to a 45 degree 
line, (i.e., an ellipse with the minor axis collapsed to zero). The "zero length" 
20 indicates that the second eigenvalue is "zero" and indicates that the eigenvector 
corresponding to the dominant eigenvalue is aligned with the major axis. The 
coherence is maximum, with a value of one. 

The situation for two traces having equal amplitudes and a 45 degree 
phase difference is shown in FIG. 7. This demonstrates how a phase shift 
25 lengthens the minor axis and therefore increases the magnitude of the second 
eigenvalue. Ihe two eigenvectors scaled by their respective eigenvalues are also 
shown. Due to the difference in these traces, the coherence is reduced to a value 
of less than one. 

Finally, in FIG. 8 both the phase and amplitude are allowed to vary. The 
30 two traces have both a 45 degree phase shift and a 2 to1 amplitude ratio. The 
resulting ellipse has a nonzero minor axis (the second eigenvalue is nonzero), 
reflecting the phase difference. In addition, the ellipse and eigenvectors are rotated 
due to the amplitude difference. Again, the amplitude and phase variations produce 
a reduction in coherence. 
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The main point of the forgoing discussion is to show heuristically that 
coherence can be expressed as a function of the eigenvalues, X, and \ 2 , and the 
eigenvectors v, and v 2 . Functionally, an expression for coherence is: 

Coherence = / (X,, X a , v,. v a ). C) 

5 This procedure for two traces can easily be extended to any desired 

number of traces. From a practical standpoint the computational burden increases 
with an increasing number of traces and limitations are imposed only by the 
available computing power. For a 3-D seismic data set. this analysis can be 
repeated over a moving or running analysis window in space and time (or space 
10 and depth), resulting in a measure of coherence at the center of the moving 
window. The output is a 3-D data set consisting of coherence values defined over 
the original data volume. 

It will become apparent from the following discussion that one advantage of 
this process is that different aspects of the data are distributed among the 
15 eigenvalues and the eigenvectors. Information, such as amplitude and phase, can 
now be analyzed and dealt with in a robust and rigorous way. Resolution, for 
example, can be improved by manipulating the eigenvalues and eigenvectors. The 
benefits can be observed visually in the computed coherence slices. Another 
important aspect of coherence values in accordance with the present invention is 
20 that they exhibit sensitivity to discontinuities in the original data and reveal subtle 
geological features, such as faults and channels. 

Returning to the process of the present invention, the first step is to obtain 
or access a set of three dimensional seismic data. Such data is in the form of 
seismic signal traces distributed over a three dimensional volume of the earth. 
25 Methods by. which such data is obtained and reduced to digital form for processing 
as 3-D seismic data are weM known to those skilled in the art Such data is 
routinely acquired by geophysical vendors who specialize in land surveys or ocean 
surveys. Such data is also sold or licensed by vendors and is generally conveyed 
or stored on magnetic tape for transfer into the memory of a seismic work station. 
30 The next step is to divide the 3-D data set into a plurality of cells or analysis 

cubes 20 (See FIG. 9). These cubes 20 perform the function of dividing or sorting 
the seismic data into groups or cells for further processing. In effect one analysis 
cube sweeps through the entire seismic data set or 3D data cube 30. Each 
analysis cube 20 comprises a stack of generally flat rectangular time layers 22. 
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For the purpose 61 simplicity, the 3-D data cube 30, the analysis cube 20, 
and the time layers 22 are shown in the form of right-angled parallelepipeds or 
cubes (generically a "ceir). Those skilled in the art will appreciate the simplicity of 
rectangular geometry in performing repetitive operations over 3-D data set. Other 
5 geometry and cell shapes are possible and may be warranted under the 
circumstances. 

Returning to FIG. 10, each time layer 22 or slab has portions of seismic 
traces t, (only one being shown to avoid cluttering the drawings) passing 
therethrough. Nine traces, laid out in a 3 by 3 evenly spaced grid are depicted in 
10 FIG. 9. Five traces in a star pattern (See FIG. 10A) or three traces (See FIG. 10B) 
may be used. A symmetric arrangement is preferred. Nine traces in each time 
layer are often better than three traces. Two traces are the least that can be used. 

To help visualize the concept of the invention, the reader is refered to 
FIG/s 11 A and 11 B. In particular, it is often difficult to visualize these concepts in 
15 nine dimensions (as is the case of 9 traces). In FIG.'s 11A and 11B, a two 
dimensional running window 24 (or analysis window) is shown with just two traces 
tt and ta contained therein. Each trace (See FIG. 1 1C) comprises a time series of 
N samples 

ti = (tn. t 12 , .... t 1N ) 

20 t^(t a ,.t» W 

The analysis window 24 of FIG. 11C is further divided into a plurality of vertically 
stacked, rectangular time layers 22. 

Before proceeding, it should be understood that in choosing window size 
and spacing, there Is a tradeoff between resolution and stability. In other words. 
25 smafl analysis windows or cubes allow higher spatial or temporal frequency in the 
resulting parameter estimate, but give less statistical stability or fewer degrees of 
freedom to those estimates. On the other hand, very large windows have poor 
resolution and tend to smear the data such that important geologic features may be 
lost 

30 Returning to FIG. 11C, the data points within each time layer 22 define 

seismic data vectors (here a 1 by 2 matrix, where N«2). As such, the two 
dimensional time layers of FIG. 11C (or the three dimensional time layer of FIG. 9) 
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form or define vector windows. Thus, nine traces (N=9) will result in a data vector 
having nine elements. 

The outer product of the data vector within each vector window or time 
layer 22 results in a N by N matrix. So a nine element data vector gives rise to a 9 

5 by 9 matrix. Adding these matricies (one matrix for each vector window) results in 
a N by N covariance matrix for the entire analysts window 24 (or the analysis cube 
20). So if the analysis cube comprises nine time layers, nine N by N matricies are 
added together to form one N by N covariance matrix. 

In accordance with the present invention, a very useful and non-obvious 

10 measure of trace coherency is obtained by computing the eigenvalues of the 
covariance matrix. In particular, the largest or dominant eigenvalue of the 
covarience matrix and the sum of the eigenvalues of the covariance matrix are 
computed. The ratio of these two numbers represents the size of the dominant 
eigenvalue relative to the sum of eigenvalues. It also indicates the variability of 

15 traces within the analysis cube. Expressed mathematically, a useful seismic 
attribute is represented by: 

(2) 

i - 1 

where K are the eigenvalues of the covariance matrix, and X, is the dominant 
eigenvalue. As such, ¥ is an indication of how well the point elements of the 
20 seismic data vectors are correlated (See FIGs. 4 and 5). 

Each and every eigenvalue of the respective covariance matrix need not be 
expressly calculated. Those skilled in the art know that methods exist for 
calculating just the dominant eigenvalue (e.g., the power method, the Rayleigh 
quotient (the fester of the two], etc). Moreover, it is also known that the sum of the 
25 diagonal elements of the covariance matrix is equal to the sum of the eigenvalues 
of that covariance matrix. For convenience, the ratio (i.e., a measure of coherency) 
of equation (2) can be assigned to the center of the analysis cube 20 or the 
analysis window 22. It should be appreciated that one advantage of using the 
dominant eigenvalue is that it tends to more directly show the variability of the 
30 traces within the analysis window. The dominant eigenvector would not be a 
measure of the variability (i.e., the coherence). 
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In performing this measure of coherency, the analysis cube 20 or analysis 
window 22 effectively sweeps (i.e.. laterally and vertically) accross the entire 3-0 
volume of FIG. 9 or all the traces of FIG. 11 A. Preferably, adjacent analysis cubes 
20 or vector windows 22 overlap one another (See windows 24, 24' and 24" of FIG. 
5 11 A). Overlap improves spatial resolution. 

The result of the analysis cube 20 sweeping over the entire 3D data volume 
30 and the assignment of the coherency measure of equation (2) is an array of 
coherency values assigned along each trace at the location of each data vector. 
As such, the 3D data volume is converted into a 3D "coherency cube." 

10 The data or coherency measures contained within the coherency cube is 

most conveniently interpreted by displaying the coherency data in the form of a 
seismic attribute map. Such a map is most often in the form of a display of those 
coherency values lying along a surface passing through the coherency cube. Two 
examples are a flat plane passing through a common horizontal time slice, and a 
15 curved surface passing through a seismic horizon line selected by a seismic 
interpreter. Another example is a line representative of geologic depositional time 
so as to capture features from the same geologic age. 

The coherency values are easily displayed for interpretation when 
displayed as a shade of gray (e.g., white indicating the highest coherency and black 
20 indicating the lowest coherency) or another color scale. (See US Patent 4,970.699 
for a -Method for Color Mapping Geophysical Data"). 

Landmark and GeoQuest interpretive workstations, for example, can be 
used to view and interpret faults and stratigraphic features by loading the 
coherency cube as a seismic volume. Such workstations are commonly used by 
25 those skilled in the art Unprocessed 3D seismic data can be conveniently loaded 
onto the workstation by means of a magnetic tape or disk which is encoded with 
instructions tor the computer to perform the above-described process. 
Visualization software (e.g., Landmarks's SeisCube software) may be employed to 
rapidly slice through the coherency cube to aid in understanding complex fault 
30 relationships. Coherency displays, including printouts in the form of seismic 
attribute maps, can reduce interpretation cycle time when used in selecting which 
seismic lines to interpret, enabling the interpreter to work around poor data areas. 
In addition, subtle stratigraphic features and complex faulting, which are not readily 
apparent on traditional seismic displays, can be rapidly identified and interpreted. 
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FIG'S. 13A, 13B "and 13C provide comparisons of the same seismic 
information displayed and processed by other processes and in accordance with 
the present invention. The differences are readily apparent. 

Of course, the process of the invention is most conveniently carried out by 
5 writing a computer program to carry out the steps just described. Such processes 
are routinely carried out on work stations previously identified. 

An elementary process flow diagram is illustrated in FIG. 12. In one 
embodiment of the invention, a software program is written in FORTRAN 77 to 
perform the process just described. 3-D seismic data 30 is read into memory. 
1 0 Based on the size and content of the available seismic data, an initialization step 32 
is performed, program parameters are initialized, data ranges are established, 
preliminary checks are made, and window sizes are set. Defaults are read or user 
selected options are read. Next processing 34 is begun. In particular, subroutines 
36 are called to sweep the data volume with an analysis cube. Within each 
15 analysis cube, a subroutine 38 computes the covariance matrix, and another 
subroutine 40 computes the dominant eigenvalues, the sum of the eigenvalues and 
the resultant coherency value. Finally, the results are combined 42 and the 
computed values are stored 44 in the form of a coherency cube. Thereafter, a 
workstation operator can access the coherency cube to display selected portions 
20 (e.g., time slice through the cube) on a CRT 46, to create a printout or seismic 
attribute map 48, to perform further analysis or to transfer to memory or to a tape 
50 for further processing elsewhere. 

Those skilled in the art are cautioned to use the method of the invention 
with due regard to what is known about the stratigraphy and geology of the region 
25 covered by the 3-D survey. 

Coherency maps have been run on several 3-D surveys. At depths of 
reasonable data quality, approximately 90% of the faults can be readily identified. 
Faults were identified on coherency maps which were very subtle on seismic 
sections, but clearly present on the coherency maps because of the robustness of 
30 the method and the map perspective of fault patterns. Since coherency maps can 
be run on uninterpreted time slices, the present invention offers a means to greatly 
accelerate mapping of the structural framework and to reveal details of fault 
relationships which would otherwise be interpreted only through tedious fault 
picking. 
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Coherence maps were generated along picked horizons 
and clearly identified salt/shale diapirs offshore 

In other locations mud and gas volcanos were clearly 
5 indicated using the process of the invention. 

Several coherence time slices showed remarkable details 
of stratigraphic features, such as abandoned river channels, mud 
flows, point bars and submarine canyons. On seismic sections, 
these features were sometimes apparent but in some cases, were 
1 0 unidentifiable even with close scrutiny. 

This invention, like that of Bahorich et al. and Marfurt et at, provides a 
method of revealing fault planes within a 3-D volume where no fault reflections have 
been recorded. Faults are often critical to the accumulation of oil. A fault may form 
a seal by cutting off a structural or stratigraphic feature so the oil is trapped against 
15 the fault On the other hand, if the fault plane contains rubble that has not been 
cemented, it may form a conduit for fluids. This may allow the hydrocarbons to drift 
up the fault plane into the feature and be trapped in it or to escape from the feature 
by drifting up the fault plane out of it Thus, fault lines can predict flow patterns in a 
reservoir and communication between injector and producing wells, for example. 

20 Seismic discontinuities can also provide the needed link to enable reservoir 

prediction between the wells and establish reservoir continuity and flow patterns 
across a field. Coherency technology can be used for finding, identifying and 
mapping of subsurface structural and sedimentological features such as faults, salt 
diapirs, unconformities, channel systems, karsting and carbonate reef fades which 

25 are commonly associated with the entrapment and storage of hydrocarbons. 
Therefore, this technology aids in the finding, extraction and production of 
hydrocarbons. In addition, it is used for identifying both shallow and deep drilling 
hazards (e.g., places where there is gas that is too close to the surface or where 
there are instabilities). Still another example is the use of the invention to search 

30 for leakage paths from known reservoirs or underground storage caverns. 

Coherency mapping of 3-D seismic is an extremely powerful and efficient 
tool for mapping both structure and stratigraphy. The new method is particularly 
sensitive to any lateral variation in wavelet character and therefore is particularly 
sensitive to the common causes of lateral variations in the wavelet (i.e., fault 
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displacement or stratigraphic variations). Thus, the subject invention encompasses 
a process, devices on which the process is recorded in the form of computer 
instructions, the product (e.g.. a map) of that process, and the manner in which 
such a product is used in the exploration of gas and oil. 

5 From the foregoing description, it will be observed that numerous 

variations, alternatives and modifications will be apparent to those skilled in the art. 
Accordingly, this description is to be construed as illustrative only and is for the 
purpose of teaching those skilled in the art the manner of carrying out the invention. 
For example, the seismic traces have been depicted as having equal spacing. 

10 Unevenly spaced traces (See FIG. 10C) may be conveniently converted into even 
spacing by interpretation. As another example, it may be useful to filter the traces 
forming the data vector to eliminate outiyers. A median filter may be used to sort 
the elements of each data vector. Edge cutoffs may be used to achieve further 
smoothing. 

15 Moreover, other algorithms may be used to measure the similarity of 

nearby regions of seismic data or to generate the "coherence cube." The 
coherency value or seismic attribute of equation (2) serves as a rather robust 
estimate or measure of signal discontinuity within geologic formations as well as 
signal discontinuities across faults and erosional unconformities. Other 
20 combinations of the eigenvalues of the covariance matrix are suggested (e.g.. 
arithmetic mean, root mean square, average, median, square root of the sum of the 
squares, square root of the product of the squares, minimum, maximum, sum. 
product etc.). Moreover, the process of the invention may be combined with other 
attributes (e.g.. AVO slope, etc.) and also applied to multi-component seismic data. 
25 Also certain features of the invention may be used independently of other 

features of the invention. For example, geologic features identified in accordance 
with the present invention can be overiayed with a velocity map to provide a means 
of cross checking the velocities. 

In addition, while coherence slice maps by themselves are very powerful 
30 mapping tools, when used in conjunction with reconnaissance mapping of 
amplitudes and dip maps, there is promise of a technological milestone in mapping 
effectiveness for the Gulf of Mexico or similar basins with 3-D seismic data readily 
available. It is believed that detailed mapping of structure and stratigraphy will be 
accelerated by mapping in a map view and less by traditional line by line picking. 
35 Interpretation in a map view of "reconnaissance" data offers significant 
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improvement in quality and quantity of interpretation. Moreover, the process of the 
invention is inherently fast Such speed aids in rapidly making bidding choices 
when concessions become available. 

Finally, it should also be understood that the principle of the invention 
5 should be equally applicable to other fields (e.g.. passive sonar, in which case the 
sensors would be acoustic and the signal sources could be hostile submarines; 
earthquake and nuclear weapon detonation detection systems, in which case the 
sensors would be seismic and the signal sources could be earthquake or explosion 
epicenters; astronomical interferometry, wherein the sensors would be 
10 radiotelescopes and the signal sources could be distant galaxies or quasars; and 
phased-array radars in which case the sensors would be the array antennae) 
wherein signals (e.g., radar, sonar, radio frequency energy, etc.) are processed to 
form images or locate changes in the structure represented by such images. 

Thus, it will be appreciated that various modifications, alternatives, 
15 variations, and changes may be made without departing from the spirit and scope 
of the invention as defined in the appended claims. It is, of course, intended to 
cover by the appended claims all such modifications involved within the scope of 
the claims. 
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CLAIMS 

We claim: 

1 . A method for the exploration of hydrocarbons, comprising the steps 

of: 

5 a) obtaining a set of seismic traces distributed over a pre- 

determined three-dimensional volume of the earth; 

b) dividing said three-dimensional volume into a plurality of 
vertically stacked and generally spaced apart horizontal time layers and arranging 
said time layers into a plurality of cells that extend laterally and vertically, each of 

10 said time layers having portions of at least two seismic traces located therein that 
define a data vector; 

c) calculating in each of said time layers of said cells the 
outer-product of said data vector 

d) combining said outer products to obtain a covariance 
1 5 matrix for each of said cells; 

e) computing in each of said cells a measure of the 
coherency of said seismic traces, wherein said measure of coherency is at least a 
function of the largest eigenvalue of said covariance matrix; and 

f) forming a seismic attribute map from a plurality of said 
20 measures of coherency of said seismic traces. 

2. The method of Claim 1, where in step (f) said map is formed by 
displaying said measures of coherency relative to a surface passing through a pre- 
determined seismic horizon. 

25 

3. The method of Claim 1 t where in step (f) said map is formed by 
displaying said measures of coherency relative to a surface passing through a pre- 
determined time line. 

30 4. The method of Claim 1, where in step (b) said cells comprise 

analysis cubes having at least five seismic traces located therein; and where in 
performing step (c) each outer product is in the form of a 5 by 5 matrix. 
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15 

5. The method of Claim 4, where in step (b) said ceils comprise 
analysis cubes that have portions of at least nine seismic traces located therein, 
and wherein said data vectors have nine elements. 

5 

6. The method of Claim 5, where in step (b) said nine seismic traces 
are arranged into a three by three grid. 

7. The method of Claim 1, where in step (b) said cells are less than 
10 100 milliseconds in thickness. 

8. The method of Claim 1 , wherein step (c) is performed in the time 

domain. 

15 9. The method of Claim 1, where in step (e) Is performed by: 

computing said largest eigenvalue of said covariance matrix, computing the sum of 
the eigenvalues of said covariance matrix, and computing the ratio of said largest 
eigenvalue to said sum of the eigenvalues of said covariance matrix. 

20 10. The method of Claim 9, wherein said sum of the eigenvalues of 

said covariance matrix is computed by forming the sum of the diagonal elements of 
said covariance matrix. 

11. The method of Claim 1, where in performing step (b) one of said 
25 two seismic traces in each ceil is located in an adjacent cell, such that said cells 
spatially overlap one another. 
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12. A method of locating subterranean features, faults, and contours, 
comprising the steps of: 

a) acquiring 3-0 seismic data covering a predetermined 
volume of the earth, said data comprising seismic traces characterized by time. 

5 position and amplitude: 

b) dividing at least a portion of said volume into at least one 
array of relatively small, adjacent overlapping, three-dimensional analysis cubes, 
wherein each of said analysis cubes contains at least three laterally separated 
seismic traces, and dividing each analysis cube into a plurality of sample intervals 

1 0 such that each sample interval defines a plurality of one by three data vectors; 

c) computing a seismic attribute for each said cube that is a 
function of the dominant eigenvalue of a covariance matrix formed from outer 
products of said data vectors; and 

d) storing said seismic attributes of said analysis cubes for 
1 5 display in the form of a two-dimensional map of subterranean features. 

13. The method of Claim 12, where in step (c) said seismic attribute is 
a function of the ratio of said dominant eigenvalue to the sum of at least two of the 
eigenvalues of said covariance matrix of said cube. 

20 

14. The method of Claim 13, where in step (c) said seismic attribute is 
a function of the ratio of said dominant eigenvalue to the sum of all the diagonal 
elements of said covariance matrix. 

25 15 . The method of Claim 14, where said seismic attribute is assigned 

to the center of its analysis cube. 

16. The method of Claim 15, wherein step (b) is performed on a 
plurality of time slices further including the step of: 
30 e) displaying said seismic attributes of successive time slices 

passing through said centers of said analysis cubes, to identify relative space and 
time invariant features. 
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17. In seismic exploration wherein 3-D seismic data comprising 
reflected seismic energy is recorded as a function of time to produce a series of 
seismic traces, and wherein a computer is used that is adapted to process such 
5 seismic traces, an article of manufacture comprising: 

a medium that is readable by a computer and that carries instructions for 
said computer to perform a process comprising the steps of: 

(a) accessing 3-D seismic data over a predetermined volume, 
said data comprising seismic signal vectors characterized by time, position and 

10 amplitude; and 

(b) ascertaining the similarity of nearby regions of said 3-D 
seismic data of said volume by: 

(1) dividing at least a portion of said data into an array 
of relatively small, adjacent, overlapping, three-dimensional analysis cubes, 

1 5 wherein each of said analysis cubes contains at least two data vectors: and 

(2) computing a seismic attribute for each cube that is 
a function of the principal eigenvalue of a covariance matrix that is formed from a 
sum of outer products of said vectors of said cube. 

20 18. The article of manufacture of claim 17, wherein said medium 

carries instructions for said computer to perform step (2) by computing the ratio of 
said principal eigenvalue to the sum of eigenvalues of said covariance matrix. 

19. The article of manufacture of claim 17, wherein said medium 
25 carries instructions for said computer to perform step (2) by computing the ratio of 

said principal eigenvalue to the sum of the diagonal elements of said covariance 
matrix. 

20. The article of manufacture of claim 19, wherein said medium 
30 carries instructions for said computer to perform step (1 ) by forming analysis cubes 

having a generally rectangular array of least five seismic traces located therein; and 
wherein said covariance matrix is at least a Ave by five matrix and is formed from at 
least three outer product matrices. 
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21. The article of manufacture of claim 20, wherein said medium 
carries instructions for said computer to assign said seismic attributes to the center 
of its analysis cube. 

5 

22. In seismic exploration wherein reflected seismic energy is recorded 
as a function of time to produce a series of seismic traces, a method comprising the 
steps of. 

(a) determining the outer-product two data vectors formed 
10 from at least two seismic traces; 

(b) forming a covariance matrix by adding said outer-products 

of step (a); 

(c) computing a seismic attribute that is a function of at least 
the principal eigenvalue of said covariance matrix of step (b); 

16 (<j) repeating steps (a) through (c) across at least a part of at 

least one time window; and 

(e) forming a map of said seismic attributes over said time 

window. 

20 23. The method of Claim 22, wherein step (c) is performed by 

computing the ratio of said principal eigenvalue to at least a partial sum of the 
eigenvalues of said covariance matrix. 

24. The method of Claim 22, wherein step (c) is performed by 
25 computing the ratio of said principal eigenvalue to at least a partial sum of the 

diagonal elements of said covariance matrix. 

25. The method of Claim 22. wherein step (d) is performed by using at 
least one seismic trace from the previous performance of step (a) and at least two 

30 new seismic traces that are located adjacent to said at least one seismic trace. 
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26. The method of Claim 22, wherein step (a) comprises the steps of 

(1) accessing 3-D seismic data over a predetermined volume 
of the earth, said 3-D seismic data seismic traces characterized by time, position 
and amplitude; and 

5 (2) dividing a portion of said volume into at least one time 

window comprising an array of relatively small, overlapping, three-dimensional 
analysis cubes that contain said at least two seismic traces. 

27. A method of seismic exploration, comprising the steps of: 

10 a) reading a 3-D seismic data set comprising seismic signal 

traces distributed over a volume of the earth; 

b) selecting at least one horizon slice from said volume and 
forming thereon cells that are arranged into laterally extending rows and columns, 
each of said cells comprising at least three seismic traces extending generally 

15 therethrough; 

c) computing for each of said ceils; 

(1) the outer-product of data vectors defined by a 
plurality of time intervals on each side of the center said cell; 

(2) a covariance matrix from said outer products of 

20 step (1); and 

(3) at least the largest eigenvalue of said covariance 

matrix; and 

(d) examining said eigenvalues of said cells along said at least 
one horizon slice. 

25 

28. The method of Claim 27, wherein step (3) is performed by 
displaying a representation of said largest eigenvalues of said cells across at least 
one horizontal time slice. 

30 29. The method of Claim 28, wherein said representation is a function 

of said largest eigenvalue of said cell and the sum of the eigenvalues of said 
covariance matrix of said cell. 
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30. In seismic exploration wherein reflected seismic energy is recorded 
as a function of time to produce a series of seismic traces, a method comprising the 
steps of: 

5 (a) arranging 3-D seismic data into relatively small, 

overlapping, three-dimensional analysis cubes that contain a plurality of seismic 
traces; 

(b) determining the outer-product of the data vectors defined 
by said analysis cubes; 
10 ( C ) forming a covariance matrix for each cube by adding said 

outer-products of step (b); 

(d) computing a seismic attribute that is a function of the ratio 
of the principal eigenvalue of each covariance matrix to the sum of all the 
eigenvalues of that covariance matrix; and 
15 ( 6 ) arranging said seismic attributes for being displayed as a 

map. 



31 . A device, comprising: 

computer readable means carrying instructions for a process comprising the 
20 steps of 

(1) reading into memory 3-D seismic data that covers a pre- 
determined volume of the earth; 

(2) digitally sorting said 3-D seismic data into an array of 
relatively small three-dimensional cells wherein each of said cells contains at least 

25 three seismic traces; 

(3) calculating in each said cells a coherency value from the 
eigenvalues of a covariance matrix formed from a plurality of outer products of said at 
least three traces; and 

(4) storing said coherency values of said cells for displaying 
30 same in the form of two-dimensional map of subterranean features represented by 

said coherency values. 
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32. The device of Claim 31, where in step (3) said coherency value is at 
least a function of the largest of said eigenvalues of said covariance matrix. 

33. The device of Claim 32, wherein said coherency value is a function of 
5 said largest eigenvalue and a sum of said eigenvalues. 

34. The device of Claim 31 wherein said computer-readable means is 
selected from the group consisting of a magnetic tape, a magnetic disk, an optical disk 
and a CD-ROM. 

10 

35. A method of locating subterranean features, faults, and contours, 
comprising the steps of: 

a) obtaining seismic data covering a pre-determined volume of 

the earth; 

15 b) dividing said volume into an array of relatively small three- 

dimensional cells wherein each of said cells is characterized by at least two seismic 
data vectors located therein; 

c) computing a covariance matrix from the outer products of 
said data vectors 

20 d) mapping a representation of the eigenvalues of said 

covariance matrix. 

36. The method of Claim 35 where in step (c) is performed by using a 
covariance matrix formed from adding a plurality of outer products. 

25 * 

37. The method of Claim 35, wherein step (d) is performed by mapping 
the ratio of the largest eigenvalue to a sum of eigenvalues. 
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38. A method of prospecting for hydrocarbon deposits, comprising the 
steps of 

a) obtaining a seismic attribute map of coherency values of 3-D seismic 
data for a predetermined three-dimensional volume of the earth, said map being 
generated by using a computer and a program for said computer that instructs said 
computer to perform the following steps: 

(1) reading said data and sorting said volume into an 
array of relatively small three-dimensional cells, wherein each of said cells has at least 
two seismic data vectors located therein; and 

(2) calculating in each of said ceils a coherency value for 
said seismic traces that is a function of the eigenvalues of a covariance matrix formed 
from outer products of said data vectors; and 

(b) using said map to identify subsurface structural and sedimentological 
features commonly associated with the entrapment and storage of hydrocarbons. 

39. The method of Claim 38. further including the step of using said map 
to identify drilling hazards. 

40. The method of Claim 39, further including the step of drilling at a 
20 location identified in step (b). 

41. The method of Claim 38, wherein said program instructs said 
computer to perform step (a)(2) by: 

(i) computing the largest eigenvalue of each covariance matrix 
25 and the sum of the eigenvalues of said covariance matrix; and 

(U) computing the ratio of said largest eigenvalue to said sum. 

42. The method of Claim 41. where in performing step (i) said program 
instructs said computer to compute said sum of eigenvalues by calculating the sum of 

30 the diagonal elements of said covariance matrix. 
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43. A seismic map prepared by a process, comprising the steps of: 

(1) accessing, by means of a computer, a dataset comprising 
seismic signal traces distributed over a pre-determined three-dimensional volume of 
the earth; 

5 (2) dividing said three-dimensional volume into a plurality of ceils 

that are arranged in time and space, each of said cells having located therein a 
plurality of data vectors; 

(3) computing in each cell a plurality of outer products formed 
from data vectors located therein; 

1 o (4) combining said outer products to form a matrix for each ceil; 

(5) computing the dominant eigenvalue of said matrix and the 
sum of the diagonal elements of said matrix; and 

(6) displaying said largest eigenvalue relative to said sum for 
each matrix of a pre-determined group of cells passing through a pre-determined 

15 surface. 



44. The seismic map of Claim 43, wherein step (6) is performed by 
obtaining the ratio of said dominant value to said sum for each said matrix of said pre- 
determined group of cells passing through said pre-determined surface. 

20 

45. The seismic map of Claim 43 where in step (2) each of said data 
vectors has at least three elements. 



46. The seismic map of Claim 43, wherein step (4) is performed by 
25 adding together said outer products. 



47. 



a) 

thereon; and 



A map for the exploration of oil and gas, comprising: 

a generally flat medium for recording visually perceptabte images 
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b) a plurality- of images on said medium that are a function of the 
dominant eigenvalue of a covariance matrix that is formed from the outer products of 
a running window of data vectors representing a 3-D seismic survey. 

5 48. The map of Claim 47 t wherein said images are a function of the ratio 

of said dominant eigenvalue to the sum of the eigenvalues of said covariance matrix. 

49. The map of Claim 47, wherein said medium is the face of a cathaode 
ray tube. 

10 

50. The map of Claim 47, wherein said images are a function of the ratio 
of said document eigenvalue to the sum of the diagonal elements of said covariance 
matrix. 

15 51. The map of Claim 47, wherein each running window comprises an 

analysis cube comprising a stack of at least three time layers; wherein each time layer 
contains a data vector therein; and wherein said data vector comprises at least three 
seismic trace elements. 

20 52. The map of Claim 51, wherein said eigenvalues are assigned to the 

center of each analysis cube. 

53. An exploration map formed by a process comprising the steps of. 

a) forming a coherency cube from data vectors of 3-D seismic data, said 
25 coherency cube comprising a three dimensional array of coherency values that are at 

least a function of the dominant eigenvalues of covariance matricies of said data 
vectors; and 

b) displaying said coherency values as an imag* on a surface in 
accordance with a predetermined transfer criteria. 



30 
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54. The map of Claim 53, wherein said coherency values are assigned to 
three dimensional coordinates that are generally coincident with elements of said data 
vectors. 

5 55. The map of Claim 53, where in step (b) said surface is a flat plane 

and said predetermined transfer criteria is that said flat plane is generally coincident 
with a time slice through said 3-D seismic data. 

56. The map of Claim 53, where in step (a) each coherency value is at 
10 least a function of said dominant eigenvalue and the sum of the eigenvalues of the 

respective covariance matrix. 

57. A device for use by a computer work station of the type used in the 
exploration of oil and gas, comprising a medium readable by a computer and 

15 containing a representation of a coherency cube, said coherency cube comprising 
measurements of the coherency of 3D seismic data, each of said measurements 
being a function of the eigenvalues of a covariance matrix formed from adding at least 
two outer products of at least two seismic data vectors. 

20 58. The device of Claim 57, wherein said data vectors are characterized 

by space and time coordinates; and wherein said measurements of coherency are 
assigned to said space and time coordinates. 

59. The device of Claim 58, wherein each of said measurements is at 
25 least a function of the dominant eigenvalue of the respective covariance matrix. 

60. The device of Claim 59, wherein each of said measurements is at 
least a function of a sum of eigenvalues. 
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A system and computer 
program for determining the location 
of a seismic event in an earth 
formation zone using seismic traces 
corresponding to the event collected 
from one or more sensors in the zone. 
A graphical data window is displayed 
of the seismic traces corresponding 
to the event, the traces representing 
observed arrival time signals from 
the sensors plotted with respect to 
time. A window of a base map of 
the positions of the sensors within the 
zone is also displayed. Displayed in 
the data window is a curve plotting 
calculated arrival time signals for the 
traces, the curve being developed from 
a user-selected one of the observed 
arrival time signals displayed in the 
data window and from a user-selected 
trial location of the event selected in 
the base map window. The program 
enables user-repositioning of the trial 
location in the base map window such 
that the displayed curve moves in 
correspondence with the repositioning 
of the trial location. The trial location 
repositioned until an acceptable 
visual coincidence of the curve is 

reached with the observed arrival time signals of the traces to thereby yield the location of the event Th<» nmemm iiw.i.„i M . • l • 
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SYSTEM FOR LOCATING SEISMIC 
EVENTS DURING EARTH FRACTURE PROPAGATION 

The present invention pertains generally to a system for detecting, recording 
and analyzing seismic events, e.g., "micro-earthquakes," caused by the propagation 
5 of hydraulic fractures occurring when fluid or slurried solid waste materials arc 
injected into an earth formation zone through an injection well, so that the fracture 
propagation may be controlled. Detection of seismic characteristics is a technical field 
in which geologists are continually in search of better tools. The invention pertains 
to apparatus and a method for quickly and accurately locating the seismic events that 
10 occur during such earth fracture propagation. 

Hydraulic fracturing of certain zones of earth formations is a now commonly 
used method to stimulate the production of hydrocarbon fluids, for example. More 
recently, it has been proposed to dispose of certain fluids or slurried solid waste 
materials into certain zones of earth formations which will contain these materials by 
15 hydraulically fracturing the zone to a predetermined extent to at least partially provide 
space for disposal of such waste material. However, a longstanding problem in 
causing hydraulic fractures is determination of the lateral, as well as the vertical, 
extent or growth of the fracture away from the point of fracture origination, typically 
a fluid injection well. In particular, the geometry of hydraulic fractures created by the 
20 injection of waste materials into a formation zone of interest is critical to avoid 
causing the injected material to flow beyond predetermined boundaries and 
contaminate aquifers or other paths which may lead the waste material into unwanted 
areas. It is also desirable to measure the growth or propagation of hydraulic fractures 
to properly control the production of fluids from formation zones of interest. 
25 An important consideration in monitoring the propagation of hydraulic fractures 

is determination of the fracture boundaries, including the radial extent of the fracture 
from an injection well on a real-time basis. This allows, for example, the injection 
of liquid and slurried solids waste materials into an earth formation to be controlled 
to prevent deposition of the materials in unwanted zones. Certain governmental 
30 authorities have jurisdiction over what earth formations and zones of certain earth 
formations may receive waste materials. The accurate determination of the 
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propagation of a hydraulic fracture on a real-time basis is thus of considerable 
importance. 

Although the detection of seismic events such as micro-earthquakes and the 
like is a relatively developed art, there has been a need to provide a system and 
5 method for measuring certain micro-earthquake events, such as those caused by the 
propagation of hydraulic fractures, to separate the actual event indicating propagation 
of a fracture from other acoustic noise in the formation, such as fluid flow generated 
noise and surface transmitted noise from machinery and the like. Accordingly, the 
process of gathering large amounts of data generated during the extension of hydraulic 

10 fractures, and the rapid analysis and display of such data in a meaningful way to 
indicate the location of an event which is part of a fracture propagation process, is of 
utmost importance to such endeavours as the disposal of slurried solids wastes and 
other fluid wastes through subterranean injection wells. 

Moreover, heretofore there has been no direct ability to monitor the growth of 

15 hydraulically induced fractures during such growth, and the extent of fracture length, 
height, width and growth rate have been assumed from pre-fracture computations for 
a predetermined set of characteristics of the earth formation in which the fracture is 
induced. These fracture models have not always been sufficiently accurate and the 
ability to actually measure the lateral as well as vertical growth of a hydraulic fracture 

20 has been a long-sought goal. An article entitled "Active and Passive Imaging of 
Hydraulic Fractures" by P. B. Willis, et al., Geophysics: The Leading Edge of 
Exploration, July 1992, describes a system for monitoring the growth of hydraulic 
fractures from one or more instrumented monitor wells wherein the seismic events are 
measured by geophones and the resultant geophone signals are later subjected to signal 

25 analysis. An article entitled "The Application of High Frequency Seismic Monitoring 
Methods for the Mapping of Grout Injections" by E. L. Majer, The International 
Journal of Rock Mechanics, Mining Science and Geomechanics, Vol. 26, Nos. 3 and 
4, pages 249-256, 1989, discusses a real-time monitoring system wherein digital data 
at sample rates greater than 50,000 samples per second are required to capture the 

30 necessary data. A sixteen channel system digitizes the data at up to 100,000 samples 
per second on each channel, saves the waveform, picks p-wave times, locates the 
events and plots the results in "close to real time using 3-D colour graphics." 
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Thus it is recognized that increased speed in accurately locating seismic events 
in three-dimensions can be achieved by interactively picking the arrival time of the 
compressional (P) and shear (S) phases of the seismic record at each detector location. 
For example, a computer program available from Los Alamos National Laboratories 
5 can be modified so that a least-squared location is determinable for each of the 
seismic events. For the large events this process is relatively rapid, taking only a few 
minutes with the aid of automated phase picking of arrival times. However, for 
smaller events with low signal-to-noise ratio, when a large number of channels are 
used, the process can take several hours and require an interpreter, such as a skilled 
10 engineer, experienced in phase identification to achieve accurate results. 

Therefore, what is needed is an improved system and method for determining 
the extent and location of hydraulic fractures propagated in earth formations, 
particularly, but not limited to, those fractures caused by the injection of fluids 
through an injection well for waste disposal. 
15 There is further a need for a system that increases the accuracy of, and reduces 

the turnaround time for, locating seismic events of small magnitude. Such a system 
may, for example, permit regulatory monitoring of fracture growth during fluid 
injection in sufficient time to impact the injection process. 

The present invention aims to provide fast turn around time and accuracy in 
20 locating the seismic events, particularly those of low signal-to-noise ratio. 

Generally, a system is provided for monitoring the extension of hydraulic 
fractures in earth formations from a fluid injection well, for example, wherein acoustic 
signals, generated and propagated through the earth formation by the seismic events 
caused by growth of the fracture, are recorded and transmitted to a processing system 
25 for display and analysis. A pattern of signals may be observed to detect the location 
of the fracture extension or propagation, and the shape of the fracture. The system 
thus permits the real-time control of fracture growth such as in the process of 
injecting fluidized wastes into certain earth formation zones. 

The location of significant events indicating the propagation of a fracture are 
30 calculated and displayed on a three-dimensional display of an injection well and 
monitor wells, if any, and the surrounding geology. The display of the seismic events 
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may indicate the accuracy of the location and predetermined limits to the accuracy of 
the location. The display may be supplemented with engineering data extracted 
through one or more computers monitoring the injection process. For example, the 
engineering data may include, but not be limited to, the fluid injection pressure into 
5 the injection well and the fluid injection rate. 

In one aspect, the invention is directed to apparatus for determining the 
location of a seismic event in an earth formation zone using seismic traces 
corresponding to the event collected from one or more sensors in the zone, the 
apparatus including a computer having a graphical user interface. The apparatus 

10 displays a graphical data window of the seismic traces corresponding to the event, the 
traces representing observed arrival time signals from the sensors plotted with respect 
to time. A window of a base map of the positions of the one or more sensors within 
the zone is also displayed. Also displayed in the graphical data window is a curve 
plotting calculated arrival time signals for the traces, the curve being developed from 

15 a user-selected one of the observed arrival time signals displayed in the graphical data 
window and from a user-selected trial location of the event selected in the base map 
window. The apparatus enables user-repositioning of the trial location in the base 
map window such that the displayed calculated arrival-time curve moves in 
correspondence with the repositioning of the trial location. The trial location is 

20 repositioned by the user until an acceptable visual coincidence of the curve is reached 
with the observed arrival time signals of the traces to thereby yield the location of the 
event. In refinements to the foregoing, mathematical algorithms to help with the 
quality of the decision, such as correlation and semblance over the data traces, may 
be implemented. The velocity model from the earth used in the calculation is stored 

25 in a database such that it may readily be modified if the data indicates a change in the 
velocity model values or details is justified. 

The apparatus also allows for refining the location of the event using a pick 
mode technique. First, the calculated arrival time curve is displayed as timing marks 
so that each timing mark corresponds to one of the traces and represents the calculated 

30 arrival time signal for the corresponding trace. The user respositions selected ones of 
the timing marks after which the apparatus performs a functional analysis on the 
calculated arrival time signals to yield a refined location of the event. The apparatus 
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may also include means for performing a polarization analysis to remove ambiguity 
in the location of the seismic event when more than one solution is possible, typically 
caused when only a limited number or distribution of earth sensors are used. 

Advantages of the present invention include providing for changing hydraulic 
5 fracture characteristics of length, height and width as well as growth rate by altering 
the pressure and rate of the injection process. In particular, the invention provides for 
monitoring the growth of the fracture height to prevent the fracture from exceeding 
a design height or from unwanted breakthrough into a particular water, oil or gas- 
bearing formation zone, for example. There are significant environmental and 
10 regulatory implications to controlling the growth of fractures which contain waste 
materials. 

Still further, the invention provides for identifying earth formation zones 
wherein significant stress contrast and strength variations in the earth formation away 
from the injection well may cause changes from a predicted extension or direction of 

15 growth of a hydraulically induced fracture. It reduces the location time from several 
hours to tens of minutes without loss of accuracy and improved ability to locate small 
events. Small events may be the result of the fracture event being far away. 
Therefore, the invention provides the potential for increased location range. 

Advantages of the event location program are that it is user friendly, and 

20 minimal training is needed for operators who are only using map location arrival time 
curve matching, as this procedure requires very little experience in phase 
identification. The trial event location solution obtained from the map location arrival 
time curve matching technique aids greatly in phase identification when subsequent 
manual picking is required. It is further useful in locating small events with low 

25 signal-to-noise ratio where identification of discrete phases is almost impossible on 
individual traces but are recognizable from a pattern on a series of multiple adjacent 
traces. 

An embodiment will now be described, by way of example, with reference to 
the accompanying drawings, in which:- 
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Figure 1 is a schematic diagram showing a waste injection well and two 
instrumented monitoring wells in accordance with the present invention; 
Figure 2 is a plan view of the diagram of Figure 1; 

Figure 3 is a detail view showing certain basic features of a geophone unit 
5 used with the geophone arrays disposed in the monitoring wells; 

Figure 4 is a schematic block diagram showing major components of the data 
recording and processing system of the present invention; 

Figure 5 is a diagram showing steps in the method of recording and analyzing 
data in accordance with the present invention; 
10 Figure 6 is a diagram of steps in the method of continuously acquiring data; 

Figure 7 is a diagram showing steps in the acquisition and recording of data 
into a shared memory and writing the data to archival tapes; 

Figure 8 is a diagram showing the steps in a method of reviewing stored data 
recorded on a tape medium; 
IS Figure 9 is a diagram showing features of a method for displaying data 

acquired and transmitted in accordance with the present invention; 

Figure 10 is a diagram showing a typical display provided for the system and 
method of the present invention; 

Figure 11 is a diagram showing another example of a display of the data 
20 recorded and analyzed by the present invention; 

Figure 12 is a diagram showing features of a method for identifying significant 
seismic events recorded by the system and method of the invention; 

Figure 13 is a diagram showing features of an interactive event timing and 
location method in conjunction with the present invention; 
25 Figure 14 is a diagram showing features of a graphical display of events 

indicating the three-dimensional location of the event with respect to a well and a 
barrier plane adjacent to a formation in which a fracture is being extended; 

Figure 15 is a flowchart showing the operation of the seismic event location 
program of Figure 1; and 
30 Figures 16-23 are representative screen displays showing the operation of the 

seismic event location program of Figure 1 corresponding to the operation described 
in Figure 15. 
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In Figures 1 and 2 the reference numeral 10 refers, in general, to a system for 
monitoring and controlling the propagation of hydraulically induced earth fractures. 
The system 10 has a data processing system 12 that includes a unique graphical 
display computer program 14 for quickly and accurately locating seismic events 
generated by the fractures. The seismic event location program 14 is described 
generally in an overview of the system 10 with reference to Figure 4 (i.e., the program 
14 designated as "FasTrak"), and is subsequently described in detail with reference to 
Figures 15-23. 

System Overview 

Referring to Figures 1-3, the system 10 includes the data processing system 
12 which is connected to a well 15 and one or more monitoring wells 16 and 18. An 
injection system 20 is connected to the well 15 for the disposal of certain liquid 
wastes and slurried particulate solids into the earth. It is understood that the injection 
well 15 extends into an earth formation zone 22 which, for example, is predetermined 
to have an in situ stress which may be less than the in situ stress of a zone 24 lying 
above the zone 22. In this way, there is some assurance that a hydraulically induced 
fracture extending from the well 15 will not break out into the zone 24 or into an 
aquifer 26, for example, lying above the zone 24. In any event, the disposal of certain 
wastes into subterranean earth formations must be controlled so that the wastes are 
contained within a predetermined formation zone and are not allowed to migrate into 
zones which might be tapped for the production of certain fluids or from which the 
wastes might otherwise migrate. 

The well 15 includes a conventional tubular casing 28 extending from a 
conventional well head 30. An open hole portion 21 of the well extends below the 
casing 28 and the casing is perforated at 32 to allow fluids to be pumped into the well 
through a tubing string 34, for example, to exit the well through the perforations 32 
and initiate a hydraulic fracture 36. The fracture 36, which is exemplary, is illustrated 
as a classic two-winged, vertically- and radially-extending hydraulically induced 
fracture. The radial propagation of "wings" 37 and 38 of the fracture 36 away from 
the well 15 is assumed to lie substantially in a plane which is normal to the minimum 
in situ compressive stress in the formation zone 22. The term "fracture" as used 
herein is not limited to the exemplary fracture 36, but may also include a type of 
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fracture wherein the principal horizontal stresses are equal or unequal. In these 
situations fracture "networks" may develop which have either a circular or somewhat 
elliptical envelope when viewed in a horizontal plane. Still further, "fractures" may 
take the form of a zone of so-called disaggregated earth material. 
5 Figure 2 illustrates the fracture wings 37 and 38 extending away from the well 

15 in opposite directions and along a line or azimuth 39 defining the assumed fracture 
extension plane. The stress field in the formation zone 22 may be such that the 
fracture extends in other directions or may, in fact, be a substantially horizontally- 
extending fracture. However, the formation zone 22 may be selected to be such that 

10 the stress field is fairly uniform and provides for the propagation of the fracture 36 
in a conventional manner as the classic two-winged, vertically- and radially- 
extending fracture. Predetermination of the orientation of the fracture plane 39 may 
be carried out using conventional stress field determination techniques such as oriented 
cores, impression packers and tilt meters, for example. Predetermination of the 

15 orientation of the assumed fracture plane defined by the line 39 is useful, although not 
mandatory, with regard to the method of the present invention. 

Each of the wells 16 and 18, as illustrated, are provided with a conventional 
metal casing 44 and a wellhead 46, although the casing may not be required in all 
instances. Each of the monitor wells 16 and 18 is also shown with an elongated 

20 tubing string 48 extending within the well from the wellhead 46 and functioning as 
a support for an array of vertically spaced-apart sensors comprising, for example, 
geophone units 50. Other arrangements are also contemplated, such as a deployment 
of sensors with the weight carried by the cable or an additional stress cable, thus 
avoiding the tubing costs. As illustrated, respective arrays 51 and 52 of geophone 

25 sensor units 50 are disposed such that they are spaced above and below the expected 
vertical boundaries of the hydraulic fracture 36. In order to enhance the acoustic 
coupling of the formation zone 22 to the geophone units 50, each wellbore of the 
wells 16 and 18 is filled with a quantity of cement 53 after the geophone arrays are 
placed in their working positions. The geophone units 50 may, for example, comprise 

30 plural 30 Hz broad band type sensors 56, such as a type SM-11, available from 
Sensor, Inc., Houston, Texas, arranged in pods or containers 55, Figure 3. Each of 
the pods 55 comprises three orthogonally-arranged sensors 56 disposed to provide one 
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vertical and two horizontal event-sensing and signal-transmitting channels. Seismic 
signals generated by the sensors 56 are transmitted via suitable conductor means 58 
to the surface and to the data processing system 12 (Figure 1). The sensors 56 may 
also comprise accelerometers or hydrophones, for example. 
5 By way of example, the injection well 15 may be characterized by 8.625 inch 

diameter surface casing and 2.875 inch diameter tubing 34 for injecting the waste 
materials into the fracture 36. A real-time bottom hole pressure gauge, not shown, 
is preferably installed in the wellbore portion 21 adjacent the perforations 32. 
Alternatively, an array of geophones may also be disposed in the injection well 15 

10 near, above or below, the perforations 32. An exemplary location is shown in Figure 
1 wherein geophone units 50a are disposed in an open hole portion 21 of the injection 
well 15 below the perforations 32. 

The monitor wells 16 and 18, by way of example, may also have casing 
diameters similar to that of the injection well. The spacing of the monitor wells 16, 

15 18 from the injection well 15 may be on the order of 80 to 150 feet and 700 to 1000 
feet, respectively. The spacing of the geophone units 50 arranged in pods 55 may be 
in the range of 20 to 30 feet along the support tubings 48. The vertical coverage of 
the geophones may be predetermined in accordance with the expected vertical growth 
of the fracture 36 and the range and wavelength of the recorded signals. Vertical 

20 coverage or extent of the geophone arrays 51 and 52 may be on the order of 750 feet, 
for example. 

Other geophone arrays may be employed, including the placement of oriented 
arrays of accelerometer or hydrophone sensors on the outer surfaces of the casings 44, 
for example. Moreover, the arrangement of the monitor wells 16 and 18 may include 

25 a third well (not shown) with an accelerometer or geophone type sensor array disposed 
therein for a more accurate determination of the propagation and geometry of the 
fracture 36 including the thickness of the fracture and the overall height. However, 
those skilled in the art will recognize that the cost of the fracture monitoring system 
10 is directly proportional to the number of monitor wells (e.g., wells 16, 18) provided 

30 and a reasonably accurate determination of the geometry of the fracture 36 may be 
obtained using one or more instrumented wells. 
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If height is a principal concern, and it usually is, a single monitor well may be 
used. Polarization analysis of the various geophone pods can lead to a unique 
solution. 

Significant seismic activity has been observed in connection with the extension 
5 of hydraulic fractures in various earth formations. The extension of hydraulic 
fractures resulting from the injection of fluids into a formation zone occurs in multiple 
relatively small extensions or "jumps," each of which causes a micro-earthquake. If 
discrete, localized micro-earthquakes occur as the fracture extends away from a 
wellbore, both laterally (radially) and vertically, the amplitude of the seismic wave 

10 generated by these activities will typically be significant enough to be observed at 
locations remote from the fracture itself. Accordingly, by sensing and recording the 
seismic activity resulting from hydraulic fracture extension at various depths in one 
or more positions with respect to the fracture extension, this seismic activity can be 
used to determine the position of the micro-earthquake and, hence, the geometry of 

15 the fracture and its location. 

In order to have a meaningful quantity of information to determine the location 
and geometry of a fracture, a substantial number of seismic signals must be 
continuously recorded from several stations over a period of time. For example, by 
recording compressional or so-called P waves, as well as shear or so-called S waves, 

20 in conjunction with arrival times, an estimate of event locations may be provided 
using known seismic or earthquake monitoring methodology including a determination 
of the compressional and shear wave velocities of the formation in question. 
Moreover, the seismic activity resulting from hydraulic fracture extension must be 
discerned from "background noise" associated with signals generated by fluid 

25 movement through the injection well 15 and within the existing portions of the 
fracture itself. 

Problems associated with the injection of waste materials into a subterranean 
formation arise if the propagation of the fracture cannot be monitored closely as this 
event occurs. In other words, if seismic activity is recorded and analyzed at a later 
30 time, the fracture may have extended out of the predetermined boundaries prescribed 
for the injection process. Accordingly, it is important to be able to closely monitor 
fracture extension or propagation with minimal time delay so that action can be taken 
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to cease or reduce fluid injection, for example, if fracture extension is indicated to be 
in a direction which might result in uncontrolled flow of the injection material into an 
unwanted zone. It is to this end, in particular, that the present invention is directed 
to provide real-time monitoring of fracture propagation which includes the treatment 
of a significant number of signals received from the geophone arrays 51, 52 described 
and illustrated herein. 

Data Processing System 

Referring also to Figure 4, signal carrying conductors 62 and 64 from each of 
the geophone arrays 51 and 52 extend from the respective wellheads 46 to the data 
processing system 12. Figure 4 illustrates components of the system 12 sufficient to 
understand the present invention when read in conjunction with the rest of the 
specification, it being understood that a description of many conventional aspects of 
the system 12 are omitted for clarity. It is understood, for example, that the system 
12 includes suitable processors and memory to implement the functions and store the 
various computer programs and data subsequently described. 

The data processing system 12 includes a suitable connector panel and channel 
selector 66 for connecting the individual conductor wires of the cabled conductors 62 
and 64 to a data recording computer 68 which includes an analog to digital converter 
70 and a central processing unit (CPU) 72. The channel selector 66 may be 
configured to permit transmission of analog signals from selected ones of the 
geophone sensors 56 of each geophone pod to the converter 70. The analog signals 
directed to the converter 70 are irregular wave signals having a variable amplitude as 
a function of time. The converter 70 is adapted to continuously acquire signals from 
the geophone arrays 51 and 52, convert these signals to digital format and transmit the 
signals to the CPU 72 for further treatment. The CPU 72 is operable to be in 
communication with a suitable recording tape disposed in a storage device on tape 
recorder 74 for receiving the continuously converted and recorded data representing 
the seismic event signals. The tape recorder 74 may be of a type manufactured by 
Metrum Information Storage of Denver, CO as their Model RSP-2150, or Parity 
Systems, Inc. of Los Gatos, CA as their Model Exabyte 8500. The CPU 72 is 
provided with suitable input/output driver programs and data storage (e.g., Exabyte 
and Metrum) to transmit the data to the tape recorders 74. 
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Alternatively, the continuously recorded data representing the geophone signals 
may be transmitted via the buffer 73 and the network 79 having a selected 
input/output driver links (e.g., Ethernet or FDDI) for transfer to a second computer 78. 
The computer 78 includes a ninety six megabyte random access memory 80 operable 
5 to receive data from the computer 68 via the network 79 (i.e., via the Ethernet or 
FDDP data link). The network 79 input/output driver programs may reside on the CPU 
72 for transmitting data to the memory 80. 

The computer 78 includes a suitable disk type data storage device 84 and 
processor capacity for storing and operating a graphical display computer program 86 
10 (designated as FracView), a seismic display and quality control program 85 designated 
Scope, an event detector program 88 designated Detect and a contingency program 90. 
It is understood that data from the contingency program 90 may be stored in a 
memory (not shown) via a network or a tape storage device 91, e.g., Exabyte. A 
second tape recorder 92 is also operable to receive suitably acquired and organized 
15 data segments from the memory 80 as will be described in further detail herein. 

The Detect program 88 is operable to identify seismic events of a particular 
characteristic and transfer the signal data from the selected event to the storage disk 
84. A third computer or central processing unit 100 (possibly multitasked with the 
computer 78) is also be provided with the system 12 and has residing thereon certain 
20 programs 14, 104, 106, 108 (described subsequently in detail) for further selecting 
signals representing seismic events of a particular characteristic and for calculating the 
location of the selected event. Moreover, the resulting analysis of selected data 
carried out by the programs 14, 104-108 may be stored on a suitable device such as 
the disk 84 for use in the FracView display program 86. Certain engineering data 111 
25 may also be transferred to the program 86 such as the rate of pumping of fluid into 
the injection well and the pressure of the fluid in the vicinity of the perforations 32. 
The computers or central processing units 68, 78 and 100 may, respectively, be 
commercially acquired from Oyo Limited of Houston, Texas, as their Model DAS1, 
Silicon Graphics, Inc., Sunnyvale, CA, as their Model Indigo 2, and Sun 
30 Microsystems, Inc. of Sunnyvale, CA as their Model SUN 2. Other possible models 
include Hewlett Packard Models V743 and HP-E1432A analog to digital converters. 
A suitable video monitor or display device 112 is also provided and operably 
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connected to the computer 78 for displaying the output of the programs 85, 86, 14 
(Scope, FracView and FasTrak). The programs 14, 104-108, 85, 86, 88 are described 
in further detail with reference to Figures 10 et seq. 

Data Acquisition and Storage 
5 Referring now to Figure 5, there is illustrated a functional block diagram of the 

major steps in acquiring and operating on the data transmitted from the monitor wells 
16 and 18. The computer or central processing unit 68 is operable to continuously 
acquire data transmitted by the multiconductor cables 62 and 64 and to prepare the 
data for transfer to the shared memory 80 and to a backup receiver in the form of the 

10 tape on the recorder 74. This process may be carried out by a continuous acquisition 
program, generally designated by the numeral 120. This program 120 may reside on 
the CPU 72. The computer 78 is also operable to read the data acquired by the CPU 
72 in segments of predetermined duration. For example, when sampling signals from 
100 channels at 2000 samples per second data segments of ten seconds duration are 

15 obtained. A unique identifier is added to each data segment and the data segment is 
then transferred to the tape recorder 92 or to the computer 78 for placement in the 
shared memory 80. The shared memory 80 is operable to read the data storage tape 
of recorder 74 or to read the data segments directly. Suitable operating programs 122 
and 124 may reside on the computer 78 for reading the data from the computer 68 or 

20 from the storage tape on recorder 74, respectively. A suitable program 126 may also 
reside on the computer 78 for transferring data from the shared memory 80 to the 
archive tapes 92. The quality control program 85 (Scope), the event detecting 
program 88 (Detect) and the seismic event location program 14 (FasTrak) are utilized 
to store data on the disk 84 for processing utilizing the display program 86 (FracView) 

25 as discussed further below. Inputs to the program 14 may also include the engineering 
data 111 and field parameters 113. 

Referring now to Figure 6, there is illustrated a functional diagram showing the 
primary steps in acquiring data from the geophone arrays 51 and 52 and preparing the 
data for transfer to the tape on recorder 74 or to the computer 78 for further treatment. 

30 The converter 70 may be operated to receive analog data signals from the cables 62 
and 64 at a selected sample rate, from a predetermined number of channels or separate 
geophone sensors, according to certain acquisition parameters and at a predetermined 
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signal length, as indicated by input block 71. For example, as many as 100 separate 
signals may be received by the converter 70 and converted to digital format at a rate 
up to 2000 samples per second, 2 bytes per word and data segments of ten seconds 
duration. The digital data, for example, may be transferred to the buffer 73, which 
5 has two buffer devices 73a and 73b, at 0.4 megabytes per second. Data is switched 
from one buffer device 73a, 73b to the other if one device is already occupied by data, 
as indicated by decision block 75. Data is transferred from the buffer 73 faster than 
the data segment length. For example, a ten second segment of data may be called 
from the buffer 73 and have an identifying header 77a or 77b placed thereon before 
10 transfer to either the tape on recorder 74 or to the computer 78 for storage in the 
shared memory 80. Transfer of data to either the tape recorder 74 or to the computer 
78 is preselected as indicated by device selection block 81, and the data transfer is 
controlled by a driver program 79 associated with the selected one of the recorder 74 
or a selected data transfer network 79 connected to the computer 78. 
15 Referring to Figure 7 (and Figure 4), the steps of reading the data generated 

by the recording computer 68 and the transferring of the data to the shared memory 
80 are illustrated. A program (not shown) resides on the computer 78 which is 
capable of accepting data from the buffer 73 by way of either of the FDDI or Ethernet 
networks. The sampling rate, word size, number of channels and record length may 
20 be further selected and transferred to the acquisition program 80 as indicated by the 
input parameters block 700. The data is stored in the shared memory 80 in sequential 
plural segments, a total of eleven, for example, during reading of the data over either 
of the networks FDDI or Ethernet. If in step 702 a network error is detected, a 
message 704 may be transferred to the display 112. If there is no error in step 702, 
25 data are transferred from the shared memory 80 to an archival tape 92 for storage and 
may be called for analysis by the programs 85, 88 (Scope or Detect) or any other 
routines 706 which are capable of being carried out by the computer 78. Data 
transferred from the shared memory 80 to the archival tape 92 is controlled by the 
steps illustrated in Figure 7 and by the Write Tape program 93. If no tape error exists 
30 and no preexisting temporary disk files exist as indicated in step 708, data are 
transferred to the tape 92. If in step 708 a tape error is sensed, an operator prompting 
signal is provided and all data arc temporarily stored on the storage disk 84 until the 
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error is cleared. If in step 708 there is no error, in step 710 a determination is made 
if preexisting temporary disk files already exist. If not the data is written to tape 92 
as shown in step 711. If in step 710 temporary disk files exist, in step 712 the new 
data are held on the disk and in step 714 written to tape before being transferred to 
the tape 92. In this way the disk is used as a temporary holding site while tape errors 
or problems are resolved by the operator, ensuring no loss of data. 

In Figure 8, it is illustrated that the computer 78 is also operable to call one 
or the other of the tape recorders 92 to read the data thereon into the shared memory 
80, as shown in step 800. A particular segment of data may be called from the 
recorders 92 by stepping either forward or back a selected number of segments, as 
indicated by the block 802 describing functions of the Read Tape program 124 (see 
also Figure 5). 

Display and Quality Control Program (Scoped 
Referring to Figure 9, the major steps and Run-Time Options 900 available 
are described for displaying data from each of the geophone signals selected. Figure 
9 describes certain features of the graphical display and quality control program 85 
(i.e., Scope). The program 85 is operable to display data segments from the shared 
memory 80 or from the tape recorder 92, save subsets of the data segments for 
redisplay and to perform certain operations on the displayed data to enhance its 
quality, ie., readability. 

For example, Figures 10, 11 and 18 illustrate typical graphical displays of the 
selected data representing signals from selected geophones and displayed by the Scope 
program 85. Figure 10 shows a display 1000 of signals from selected geophones at 
selected positions 1-8 in the arrays 51 or 52 versus time. The geophone providing 
a signal of significant magnitude at the earliest time from the onset of the 
measurement period is indicative of the position of an event which typically represents 
the extension of the fracture 36. For example, if geophone position number 4 in 
Figure 10, at a predetermined position in the array 51 or 52, receives a signal 
indicating a micro-earthquake before any of the other geophones in the array, this 
could be correlated with the height of the fracture extension event. Since the position 
of each of the geophones in the arrays 51 and 52 is known, the geophone receiving 
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a particular seismic signal earlier than any other geophone would indicate the 
approximate depth of the event, such as a fracture propagation event. 

Accordingly, data may be processed by the steps of Figure 5 and displayed by 
the Scope program 85 very rapidly with no further analysis to provide real time 
5 display of the location (depth) of the outermost edge of a fracture. In this way also, 
the geophone arrays 51 and 52 are used to sense not only the vertical growth of the 
fracture but also its lateral or radial growth from the injection well. In regard to the 
latter determination, since the location of the geophone arrays 51 and 52 is known, 
and the acoustic velocity of both compressional wave propagation and shear wave 
10 propagation can be determined for a particular formation, the position of the fracture 
growth event from the injection well can also be determined. Figure 1 illustrates how 
geophones 50b and 50c in the respective monitor wells 16 and 18 sense the micro- 
earthquake associated with the extension of the fracture 36 at the point 36a before any 
of the other geophones in the arrays 51 and 52 receive the seismic signal. 
15 Figure 11 illustrates another display 1100 that may be obtained with the 

transmission and manipulation of the data obtained from the shared memory 80 and 
depicted on the monitor or display 112 by the program 85. The display of Figure 11 
indicates the signals received by the geophones of a particular pod 55 where the 
geophones or channels are identified by numerals 1, 2 and 3. The irregular wave-type 
20 signal illustrated on geophone 2 is of significantly greater amplitude than that received 
on geophones numbers 1 and 3 from the same event. Knowing the orientation of 
geophone number 2 with respect to the injection well 15 may be useful in indicating 
the direction and distance to the event, for example, relative to the location of the 
geophone location of corresponding traces 1, 2, 3. The display of Figure 11 is also 
25 a display of signal amplitude versus time with the time scale indicated covering 
approximately 0.50 seconds. 

Data operated on by the computer 78 using the Scope program 85 may be 
displayed in real time as acquired by the memory or retrieved from a tape on the 
recorder 92. In any case, as a data segment is replaced by a succeeding data segment, 
30 the signals displayed, such as indicated by the displays of either Figures 10, 11 or 18, 
will be replaced by the displays of succeeding signals. These signals, while displayed, 
may have their display modified to amplify both the ordinate and abscissa values, for 
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example. The signals displayed on the display 112 under the operation of Scope 
program 85 may be saved in any subset desired, the display may be reset, and the 
signal displayed may be maintained on the display 112 by an interrupt-'Tause" 
function of the Scope program 85. Still further, a mouse selectable zoom window 
5 may be provided to capture and amplify a portion of the data as displayed in Figure 
11, for example. Referring to Figure 11, there is illustrated a window 85a for part of 
the signal trace 85b wherein the portion of the signal trace within the window 85a 
may be displayed by itself and subjected to the features of the Scope program 85, that 
is the amplification of the ordinate and abscissa. Still further, the Scope program 85 
10 is operable to automatically adjust the time annotation indicated in Figure 11, for 
example, provide for labelling individual traces and for highlighting the channels 
which have received a signal or compressed display signals for a particular data 
segment. 

Event Detector Program fDetecfl 

15 Figure 12 illustrates some of the features of a micro-earthquake or seismic 

event detection method utilized in the present invention. The routine carried out by 
the event detector (i.e., Detect) program 88 is based on algorithms described in a 
treatise entitled "ASP: An Automated Seismic Processor for Micro-Earthquake 
Networks", McEvilly and Majer, The Bulletin of the Seismological Society of 

20 America, Volume 72, No. 1, pp. 303-325, February, 1982. The Detect program 88 
is operable, for selected channels or signals to identify valid events based on the time 
and amplitude of the signals received by the geophone sensors of each pod and each 
array. The parameters, which are selectable to "screen out" invalid events and to save 
so-called valid events, are identified in Figure 12 and are further explained in the 

25 above reference. Each memory segment may be analyzed as to how many, if any, 
valid events are sensed and those valid events which are identified are then given a 
unique identifying header and forwarded to a storage medium such as the disk 84. As 
indicated in Figure 2, the Detect program 88 may also receive data from a source 
other than the shared memory 80. 

30 Figure 13 depicts an alternative event location program (tPick) that utilizes an 

error minimizing method of determining the position of an event based on seismic 
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arrival times at the monitor wells 16 and 18. Still further, the seismic event may be 
located using an algorithm developed by the U.S. Geological Survey. 

Figure 14 shows one method of displaying the located events using one or 
more of the location methods described above and in detail below. The display 112 

5 shows, in Figure 14, the injection well 15, the monitor wells 16 and 18 and the upper 
boundary layer 24a defining the zone 24. Event locations identified by the Detect or 
Scope programs 88, 85 are located and displayed at positions 113, which may take the 
form of small ellipsoids using the calculation methods described. The event location 
program 14 (FasTrak) is operable to select the number of signals identified by the 

10 Detect or Scope programs 88, 85 for treatment by the event location calculation 
method selected. 

The display program 86 (FracView) is operable to provide the three- 
dimensional display indicated in Figure 14. Well trajectories, the earth's surface and 
downhole geophone locations, as well as the located seismic events indicated by the 
15 numerals 113 may be displayed in a three-dimensional display that can be zoomed, 
translated or rotated in any direction and with or without directional illumination. The 
seismic events displayed at 113 have a size scaled to the uncertainty of location, that 
is small ellipsoids indicating greater confidence in the hypo central location. The 
display is particularly effective in communicating the complicated geometry of 
20 fractures extending from a disposal or injection well to persons unfamiliar with data 
related to seismic events. 

Event Location Program (FasTrak') 
Figures 15-23 illustrate the features of the seismic event location (i.e., 
FasTrak) program 14. The program 14 is utilized to determine the location of seismic 
25 events corresponding to an extension of a fracture within an earth formation zone 22 
of interest (Figure 1). The program 14 significantly reduces the time required for 
locating seismic events, particularly smaller events with low signal-to-noise ratio 
using less-skilled staff. As will be described, the program 14 eliminates the necessity 
to manually pick individual phase arrivals, resulting in greatly reduced analysis time. 
30 The program 14 is a refinement of the program tPick (not shown) previously described 
in the parent of the present specification. 



WO 98/09182 



PCT/GB97/02345 



- 19 - 

Figure 15 is a flowchart of the steps followed in the program 14 flow. 
Figures 16-23 illustrate the graphical user interface of the program 14 as presented 
on the display 112. It is understood that the program 14 is stored and executed on the 
processing unit 78 but may alternatively be stored and executed on any other 
computer. Suitable programming instructions (not shown) control the operation of the 
program 14 as described herein. Since the instructions used to implement the 
functions of the program 14 may be understood and written by one skilled in the art 
in possession of the present disclosure, they are not described in further detail herein. 

As shown in Figures 15 and 16, execution the program 14 begins in step 1500 
by the display of a data window 1602 in a program screen 1600. The data window 
1602 is a display of seismic traces plotted vertically in increasing depth order. The 
seismic traces arc plotted showing signal amplitude (on the vertical axis of each trace) 
versus time (on the horizontal axis). Each of the traces has a reference number as 
shown and the order of the traces displayed may be chosen from a lookup table. 

The seismic event data corresponding to the displayed data traces originate 
from the earth formation zone 22 of interest and correspond to an extension of a 
fracture within the zone. As described previously, data signals are received by the 
processing system 12 from the sensors 56 (Figures 1 and 3), with the sensors being 
disposed in a predetermined array within the formation zone 22. The data signals are 
transmitted from the sensors 56 at a selected sampling rate to the A/D converter 70 
within the recording computer 68 to produce corresponding digital data signals (Figure 
4). The digital data signals are stored in the recorder 92 and the memory 80. 
Typically, the data signals are organized in data segments containing a preselected 
number of data signals for a predetermined time period, e.g., ten seconds. 

Step 1500 represents input of the seismic event data to the program 14. Before 
being input to the program 14, the seismic event data is optionally (and preferably) 
screened using the Detect program 88, described previously. The data segments are 
screened (by conditional parameters typical of those described by the Detect program 
88) to select only those data segments which have a predetermined amplitude for a 
predetermined time (thus representing an event of interest). This is done because, as 
described previously, the system 10 is continuously monitoring the earth formation 
zone 22, thus producing large volumes of data. The Detect program 88 runs in the 
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background to automatically track noise levels of the signals such that when a seismic 
event of interest occurs, i.e., when there is an anomalous event, a recording is made 
(including a sufficient time window prior to, and after the event or "tail") for storage 
and review by the program 14. In the recording of data using the Detect program 88, 
it is understood that since data is being collected from the system 10 in a stream of 
about 0.4 megabits per second on a high speed data channel (conductors 62, 64 
through channel selector 66, Figure 1), it is easy to miss events of interest and much 
unnecessary data must also be stored. Hence compression of data occurs wherein the 
Detect program 88 operates as a pattern recognition tool to store only events of 
interest for review by the program 14. It is understood, however, that all of the data 
is stored on tape for later interpretation. Note that the Scope program 85 may be used 
to identify events of interest and the output is stored on the disk 84. 

Upon receiving the seismic data from 'the earth sensors 56 as filtered using the 
detect program 88, the event location program 14 visually displays the data segments 
composed of the seismic traces and performs quick and accurate processing, as 
described in detail below, to determine the location of the seismic event 
(corresponding to the location of an extension of the fracture) in a minimum time after 

occurrence of the event. 

Still referring to Figures 15 and 16, the data window 1602 plots an illustrative 
example of data traces plotted as a function of time, the traces corresponding to a 
seismic event of interest to be located using the program 14. Each active trace 
includes a first, smaller-amplitude recording of a compressional or so-called P wave 
followed by a larger-amplitude recording of a shear or so-called S wave. Additional 
traces of auxiliary date (along the bottom of Figure 16) may be displayed in the data 
window 1602. In this case a large amplitude signal from a strong motion detector are 
displayed, the data for which may or may not be used in the event location procedure 
described below. 

In step 1504 the format of the data window 1602 is enhanced by the user to 
better view the seismic traces using options described in the Scope program 85. A 
display window 1606 appearing on the screen 1600 to the left of the data window 
1602 operates as a control box that presents a selection of program options for 
changing the for- it. 
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The options may be selected by entering computer keyboard commands (using 
alphanumeric or function keys, as indicated) or by using a pointing device (e.g., a 
mouse or trackball instruction pointer) to make the selection. The format options 
include: increasing the spaces between the traces; changing the time scale; 
increasing/decreasing trace gain; plot zooming; windowing portions of the data 
display; and saving a window to disk or to a printer. In the display window 1606, 
representative command options are shown. "Open" is used to present an open 
display where traces at a common location (x,y,z) are each plotted as separate traces, 
while "Closed" is used to present a display where traces from a common point 
location are each plotted on the same trace line but in different colours. Function keys 
F1-F4 are used to increase or decrease the X and Y scales of the traces shown in the 
data window 1602. Function key F8 returns the display in the data window 1602 to 
a default display. F9, F10 increase and decrease, respectively, the gain of the signal 
amplitudes of the traces. Fll, F12 increase and decrease, respectively, the zoom, i.e., 
the magnification, of the traces in the data window 1602. The colour Khaki is used 
for clarity to indicate which command options require coincidental use of the shift key 
to complete the command selection. 

In step 1506 a base map showing the positions of the detectors, i.e., the sensors 
56, is displayed on the screen 1600 in a map window 1604. The base map is a plan 
(aerial) view of the locations of the receiver sensors 56 relative to the injection well 
15. The map window 1604 may be sized, moved and presented in overlapping 
window formats as desired by the user and as necessary to suit the workstation 
environment. 

In step 1507 a first estimate of the event depth is made by selecting and 
highlighting a trace in the data window 1602 (Fig. 17, discussed below). A readout 
of the depth is then available in the map window 1604. 

In step 1508 the seismic event is located using a "map-mode" technique. The 
map-mode technique uses curve matching to achieve a visual agreement between a 
curve plotting the arrival times of the traces (P waves or S-waves or both) calculated 
from a trial location, and the actual recorded data of the seismic traces shown in the 
data window 1602. The visual agreement is achieved by interactively moving the trial 
location curve so that the calculated arrival times represented by the curve visually 
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agree with (i.e., they coincide to a best approximation in the opinion of the user) the 
displayed arrival times of the recorded seismic traces in the data window 1602. 

Referring also to Figures 17 and 18, the steps in performing the event location 
using the map mode technique of step 1508 are illustrated. A first estimate of the 
5 event depth is made by selecting and highlighting a trace in the data window 1602. 
In Figure 17 the trace "58 59 60 M is highlighted (typically by colour contrast) to 
identify the trace in the display window 1606 by clicking the mouse pointer over one 
of the trace annotations, say, for example, "58" whereupon the trace 58 will be 
highlighted in bright pink. 

10 In this case the compressed (closed) mode display option from the display 

window 1606 is being used to display three traces (58, 59, 60) at similar ground 
locations, such as generated from a three-component geophone pod (Figure 3), all 
being displayed on the same trace axis in three different colours in the data window 
1602. The selected trace is preferably (but not necessarily) the trace with the earliest 

15 phase (compression wave) arrival time because the depth of this trace is a good first 
estimate of the source depth of the event. A timing mark 1702 is placed on this trace 
at a control point as shown to establish a relative time for this event. While in the 
present instance the timing mark 1702 and its corresponding control point are selected 
for a compressional wave arrival time, the selection may be either for the 

20 compressional or the shear wave arrival time (the left mouse button being used for 
selecting the P wave, the right mouse button being used for selecting the S wave). 

Next, the cursor is moved to the base map window 1604 and a trial location 
1802 is selected and plotted by the user for the seismic event as a red dot (Figure 18). 
The numerical representation of the trial location 1802 is printed in a lower portion 

25 of the base map window 1604. With the depth (Z) established by trace selection using 
the timing mark 1702 control point, and the X, Y position of the trial location 1802 
provided in the base window 1604, it is possible to calculate arrival times at all other 
geophones, i.e., for all other traces. This is done by the program 14 and a calculated 
arrival time curve 1804 is drawn from the originally selected trace (the so-called 

30 "seeded" trace) to all other geophone locations, i.e., to all other traces in the data 
window 1602. 
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As seen in Figure 18, the calculated arrival time curve 1804 generated from 
the selected trial location 1802 and the timing mark 1702 control point does not 
coincide exactly with the recorded, i.e., observed, compressional and shear wave 
arrivals for the seismic traces displayed in the data window 1602. Accordingly, the 
user interactively repositions (using the mouse in a "click and drag" fashion, for 
example) the trial location 1802 in the map window 1604. The X-Y direction of the 
trial location 1802 is repositioned with the mouse and the depth (Z direction) of the 
trial location is repositioned with the up and down arrow keys, for example. The 
display of the calculated arrival time curve 1804 in the window 1602 is modified in 
real time during the interactive repositioning of the trial location 1802 such that the 
curve 1804 is moved as the trial location is moved. Thus the user can interactively 
reposition the curve 1804 by moving the trial location 1802, so that the curve visually 
coincides with the observed compressional and shear arrival times of the traces. Even 
greater accuracy is achieved by changing the data window 1602 and the base map 
window 1604 scales (e.g., by zooming in). 

Once an acceptable match between calculated and observed data is achieved 
by coincidence of the curve 1804 with the observed arrival times of the traces in the 
window 1602, the accepted location 1802 with its calculated arrival times at each data 
channel represented by the curve 1804 are then stored as a solution in a database as 
indicated in step 1514. This assumes that no greater accuracy in the event location 
solution is required as indicated in step 1512. The stored solution may then be 
displayed in three-dimension using the FracView visualization program (Figure 14) 
or otherwise processed. The location is judged for its acceptability and origin. If it 
is attributable to noises generated by normal oil field operations the solution is 
discounted. If not, the solution is assumed to be created by fracturing operations. A 
decision is made whether the fracture is growing as expected in an acceptable manner; 
if not, a decision is made to cease fracture generation. 

The process of event location using the map mode technique of step 1508 by 
selecting traces and moving the trial location 1802 to fit the curve 1804 to the 
observed arrival times of the traces (hereinafter referred to as the "map location 
process") is extremely rapid and can be accomplished with reasonable accuracy (e.g., 
within less than fifty feet) with minimal training of the operator. 
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Referring to Figures 15 and 19, further refinement in the location of the event 
is achieved with the program 14. In step 1510 a determination is made that greater 
accuracy in event location is required by the user than can be achieved with the map 
location process just described. In some cases a more accurate location of the seismic 
5 event can be obtained by modifying the shear and compressional arrival times at 
individual recording (trace) locations, because some sensor locations have varying 
delay times due to local conditions such as cement thickness, coupling or the like. In 
addition, due to changes in the rock formation due to the presence of the fracture, 
more complicated or changed velocity may be required. These can be introduced here. 

10 To compensate for these anomalies, in step 1512 the program 14 executes a 

pick mode technique to refine the location of the event by picking and adjusting the 
calculated arrival times of individual traces determined from the curve 1804 in the 
previously described map location process. Once the arrival times of individual traces 
are adjusted, a functional analysis such as a least-squares fit analysis is performed on 

15 the arrival time data to produce a more accurate event location solution (i.e., a solution 
that better matches the data). Other well-known event location algorithms (other than 
a least-squares fit analysis) may likewise be used on the adjusted arrival time data to 
determine the event location, as will be recognized by those skilled in the art. 

As shown in Figure 19, the pick mode is initiated by entering an instruction 

20 in the Pick Locate instruction box 1902. The calculated arrival times, determined 
from the map location process and represented by the moved curve 1804, are 
converted to discrete picked timing marks 1904 for each of the traces. The timing 
marks 1904 (representing the previously calculated arrival times), or selected ones of 
them, are then manually moved by the user to locations that more accurately match 

25 the observed arrival time data. 

Once the moved timing marks 1904 are in the locations desired by the user, 
an instruction is entered by the user to use the arrival times thus selected for each 
trace (by the moved timing marks 1904 and possibly others left unchanged) to locate 
the seismic event position and occurrence time. More specifically, an instruction in 

30 the Pick Locate box 1902 is entered to accept the picks, whereupon the program 14 
then performs a functional analysis (such as a least squares fit analysis) to determine 
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the new and refined event location based upon the adjustments to the arrival times of 
the traces. The refined event location is shown as described below with respect to 
Figure 20 and is also displayed in the map window 1604 as a spot of white colour. 

As illustrated in Figure 20, once the functional analysis (i.e., the inversion) is 
performed to fit the data to the newly picked arrival times, the data window 1602 
displays the picked timing marks 1904 (in green, for example) and also displays newly 
calculated arrival timing marks 2002 (in black bold dashes, for example) 
corresponding the least squares fit (or other functional analysis) performed on the trace 
data. The picked timing marks 1904 and the newly calculated arrival time marks 2002 
are presented in differing colours on each trace for comparison purposes. The user 
may then make a judgment on the acceptability of the agreement between the newly 
calculated arrival times (marks 2002) and the observed and picked arrival times 
(marks 1904). Frequently, minor adjustments to the picked arrival timing marks 1904 
are necessary, and the process is repeated accordingly, until some acceptance criteria 
is passed. The criteria is, for example, when the root-mean-square (RMS) of the 
residual of the mismatch at each trace is below a set value (such as 6 ms or 50 feet). 

As illustrated in Figure 21, a tabulation of the statistical analysis of the 
agreement between the picked and newly calculated arrival times may be generated. 
It may be printed or displayed in a window on the screen 1600 in real time during the 
inversion process and may be stored for further evaluation. The Figure 21 table, 
indicated by the reference numeral 2102, lists the P-Arrival, S-Arrival and P and S 
Computed for each Station/Channel (P/S). The Station/Channel indicates the station 
(i.e., the physical ground position, e.g., 1001), the receiver (the receiver upon which 
trace is recorded, i.e., rid=l, and the channel (or trace) number, i.e., tid=l). The P- 
Arrival represents the P wave (compressional) wave arrival time (e.g., 0.253), followed 
by the difference (in parenthesis, e.g., -0.882) between the newly calculated P wave 
arrival time (as determined by the functional analysis, e.g., the least squares analysis 
represented by the black timing marks 2002) and the picked P wave arrival time 
(represented by the picked and moved timing mark 1904). The S-Arrival represents 
the S wave (shear) wave arrival time (e.g., 0.352), followed by the difference (in 
parenthesis, e.g., -2.150) between the newly calculated S wave arrival time (as 
determined by the functional analysis, e.g., the least squares analysis represented by 
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the black timing marks 2002) and the picked S wave arrival time (represented by the 
picked and moved timing mark 1904). The P/S Computed column represents the 
newly calculated arrival times (P wave and S wave) as calculated by the functional 
analysis or the event location determined by the inversion, as represented by the 
5 calculated timing marks 2002. At the bottom of the table 2102 is an indication of the 
quality of the agreement between the picked and newly calculated arrival times, which 
in the present case is a "5," indicating poor quality (i.e., the differences between newly 
calculated and picked arrival times are large). The process can then be repeated of 
picking and performing the functional analysis until a good agreement is achieved. 

10 Various statistical analyses may be diplayed here after a solution of acceptable 

quality (Q=l) is determined. These include RMS agreement calculations and 
ellipsoids of uncertainty for each event. 

The result is an improved location of the seismic event using more detailed 
analysis of the arrival times, a possibly more sophisticated model of the velocity, and 

15 a more quantified analysis of the quality of the solution determined. 

If there is no ambiguity in the event location (discussed further below), as 
determined in step 1516, the event location solution as finally determined is stored in 
step 1514 in a database for further analysis. 

Figure 22 illustrates the performance of a polarization (hodogram) analysis 

20 on the arrival time data. In some receiver (geophone sensor) geometries there exists 
ambiguity in the calculated event location based upon arrival time alone, because 
equally acceptable event locations may be determined using the foregoing methods. 

As shown in step 1518 (Figure 15), a polarization analysis is used to remove 
the ambiguity in the location of the event. First, data collected to determine the 

25 velocity of the formation medium in the zone 22 is used to determine the orientation 
of the geophone sensors 56. The sensors 56 are typically three orthogonal sensors 
mounted close together in a pod 50. The orientation of each sensor 56 is stored in a 
database along with the location of each sensor and used as an input parameter file 
for the program 14. Optionally, data for certain traces may be modified wherein the 

30 arrival time signals for the three sensor components arc mathematically rotated into 
a principal coordinate system, usually north, east and vertical depth. Stated another 
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way, a mathematical vector multiplication is performed on the sensor signal data to 
indicate what would be the equivalent arrival time signal from each sensor component 
if the sensor component was pointed in the directions of the principal coordinate 
system. This is done for consistency of display purposes between the different sensor 
pods using a known calibration position (e.g., such as from a point of an explosion). 

Referring to the data window 1602 (Figure 22), a trace (e.g., 38, 39, 40) is 
selected and highlighted. A timing mark is indicated for the trace. Next, a time 
window is defined around the selected timing mark, the size and shape (symmetric, 
asymmetric) being retrieved from a parameter file. 

The program 14 then computes the polarization, i.e., the direction of arrival for 
each of the arrival time samples within the time window 2204. It is understood that 
different arrival time samples will be received from each of the different directional 
sensor components (traces 38, 39, 40) and each such trace of arrival time samples will 
be highlighted in a different colour. 

Next, vectors representing the directions of the arrival time samples are plotted 
in the map window 1602 as indicated by vectors 2202. The direction of the vectors 
2202 is determined using the arrival time sample and the orientation data stored in the 
database indicating the direction that the sensor is pointed in. Each of the vectors 
2202 represent what orientation the seismic source came from, as indicated by the 
directional information of the geophone sensors and the amount of energy (amplitude) 
in each of the traces. The vectors 2202 are thereby drawn that indicates the direction 
of the arrival time signal based on the amplitude of the sample and the direction (or 
orientation) of the geophone sensor. The length of each vector 2202 is drawn 
proportional to the amplitude of the arrival amplitudes at each time sample. If, for 
example, there was a purely polarized arrival all the orientations derived from samples 
within the data window would point in the same direction. 

The direction and length of the vectors 2202 thus represent the motion of the 
ground at the geophone. The motion is produced by the passing seismic wave 
initiated by the rupture event. The particle motion can be visualized on the map 
window by plotting the vectors for each sample of the recorded signals from the 
geophone. Data from all three components at a single location are necessary to 
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calculate this information and makes use of the geophone orientation data stored in 
the input files. 

The number of time samples, and hence vectors 2202, depends on the size of 
the window 2204 and typically includes about 5-10 time samples. Thus the vectors 
5 2202 are lines representing a direction defined by the amplitude of the arrival time 
signal in each of the directions defined by the coordinate system, the length of such 
line being proportional to the amplitudes of the arrival time signals. 

Once the polarization direction vectors 2202 are drawn, removal of ambiguity 
in the event location solution is accomplished by comparison of the vector 2202 with 
10 the trial location 1802 (Figure 18). The trial location 1802 is able to be then moved 
along a path defined by the lines, such that ambiguity in the event location is 
removed. When the trial location 1802 is moved to be on the polarization vector 
2204, the event location is considered to be accurately determined with the least 
ambiguity. 

15 As the trial location is moved to be on the vector 2202 the curve 1804 changes 

and an optimal solution is achieved when the calculated and observed arrival times 
coincide. Several traces are usually examined for consistency in polarization. 

Referring again to Figure 15, after performance of the polarization analysis 
(unless such is not needed in step 1516), the event location solution as finally 

20 determined is stored in step 1514 in a database for further analysis In step 1520, the 
event location solution is passed to the Display program (FracView) 86 (Figure 14), 
for viewing in three-dimension to determine if a problem exists or not. For example, 
seeing the event location solution in three-dimensional space (Figure 14) enables the 
user to determine, e.g., if the event (i.e., a fracture) is above the control horizon 24a 

25 or near a well or if some other condition exists that would result in an environmental 
hazard. In step 1522 a decision by the user based on the review is made and if the 
three-dimensional review indicates that corrective action is not required, then for that 
event the procedure ends and further monitoring of other subsequent events can occur. 
However, if in step 1522 corrective action is indicated, then in step 1514 the pumps, 

30 i.e., the injection system 20, is shut down so that the integrity of the zone 24 is 
maintained. 
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Figure 23 illustrates an example of a data window 1602 for traces which have 
very weak arrival time signals in a high noise background from an actual survey. 
Using the program 14 on the survey permits extremely weak events to be located that 
would not otherwise be easily accomplished except by a highly trained engineer. 
5 Thus the system 10 and particularly the program 14 offers several advantages. 

It reduces the time to determine the location from several hours to several minutes 
without loss of accuracy. The program 14 is user friendly, and minimal training is 
needed for operators who are only using map location arrival time curve matching, as 
this procedure requires very little experience in phase identification. The trial event 
10 location solution obtained from the map location arrival time curve matching 
technique aids greatly in phase identification when subsequent manual picking is 
required. Thus the program 14 continues to offer options using operator picked arrival 
times with greater accuracy at a higher speed than routinely done due to the map 
location trial solution. 

15 Software for the completely automated picking of arrival phases may also be 

utilized with the program 14. This would further decrease the time needed to derive 
a map location for a given event under certain signal-to-noise conditions. Further, 
deriving locations from sensors in a single monitor location is possible. It also has 
been demonstrated that a single 3-component sensor can provide a unique location 

20 using P and S arrival times together with a polarization analysis of the type described 
above. It is understood that it is not possible to locate events using a least square 
event location code from a single vertical well because of non-uniqueness. Increasing 
the number of monitor wells required in such an inversion negatively impacts the cost 
of a project, and therefore polarization analysis, combined with the rapid calculation 
25 of arrival times, allows for a single monitor well location of events at the least field 
cost. 

Those skilled in the art will recognize from the foregoing description that a 
unique system and program for seismic event location corresponding to a hydraulic 
fracture extension or propagation in a predetermined earth formation zone of interest 
30 may be carried out with the invention as described above. Although a preferred 
embodiment has been described in some detail, those skilled in the art will also 
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recognize that various substitutions and modifications may be made without departing 
from the scope and spirit of the invention. 

The preferred embodiment may be implemented as a computer workstation 
with appropriate software. It will be appreciated that certain functions may be 
5 implemented in hardware. 
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1. A computer system for determining the location of a seismic event in 
an earth formation zone using seismic traces corresponding to the event collected from 
one or more sensors in the zone, the computer having a graphical user interface, the 

5 apparatus comprising: 

means for displaying a graphical data window of the seismic traces 
corresponding to the event, the traces representing observed arrival time signals from 
the sensors plotted with respect to time; 

means for displaying a window of a base map of the positions of the one or 
10 more sensors within the zone; and 

means for map-locating the event by displaying in the graphical data window 
a curve plotting calculated arrival time signals for the traces, the curve being 
developed from a user-selected trial location of the event selected in the base map 
window, including means for repositioning the curve to coincide visually with the 
15 observed arrival time signals of the traces and thereby indicate the location of the 
event based on the changes in the calculated arrival time signals caused by the 
repositioning of the curve. 

2. The apparatus of claim 1, further comprising means for repositioning 
the curve wherein the curve is graphically moved in correspondence with movement 

20 of the trial location in the map window. 

3. The apparatus of claim 1 or 2, further comprising means for performing 
the repositioning with an instruction pointer operated by the user. 

4. The apparatus of claim 1, 2 or 3, further including means for pick- 
locating the event to further refine its location after the map-locating, comprising: 

25 means for displaying timing marks at the positions of the curve on the traces; 
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means for repositioning the timing marks at picked ones of the traces, the 
timing marks being repositioned to coincide more accurately visually with the 
observed arrival time signals, and 

means for performing a functional analysis on the calculated arrival time 
5 signals following repositioning of the timing marks to yield a refined location of the 
event based on the changes in the calculated arrival time signals caused by the 
repositioning of the timing marks. 

5. The apparatus of claim 4, wherein the means for performing a 
functional analysis is arranged to perform a least-squares statistical analysis. 

10 6. The apparatus of claim 4 or 5, further including means for automatically 

picking the traces for which the timing marks are to be repositioned. 

7. The apparatus of any preceding claim, further including means for 
providing options in formatting of the display of the seismic traces in the graphical 
data window. 

15 8. The apparatus of claim 7, wherein the formatting options include one 

or more of plot zooming, changing the trace gain, changing the time scale, and 
changing the spacing between traces. 

9. The apparatus of claim 1, wherein the map locating means further 
comprise: 

20 means for selecting a trace in the graphical data window wherein the trace 

becomes highlighted; 

means for selecting an observed arrival time of a trace to provide a timing 
mark through which an analytically predicted time of the curve passes through, such 
that all calculated arrival times of the curve are determined relative to the timing 

25 mark; 

means for selecting a trial location of the seismic event in the map window 
wherein the selected trial location is highlighted; 
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means for displaying the coordinates of the location of the event as determined 
by the curve repositioning; and 

means for storing the calculated arrival times and the event location as 
determined by the curve repositioning. 

10. The apparatus of claim 4, wherein the pick-locating means further 
comprise: 

means for converting the curve to timing marks representing the portions of 
the curve corresponding to the calculated arrival time signals; and 

means for highlighting a trace for which timing mark repositioning is to occur. 

11. The apparatus of claim 10, wherein the pick-locating means further 
comprise: 

means for generating a tabulation of the extent of agreement between the 
repositioned calculated arrival times and the observed arrival times as indicated by the 
functional analysis performed. 

12. The apparatus of any preceding claim, further including a polarization 
analyser to remove ambiguity in the location of the seismic event, the one or more 
sensors for collecting the event data including differently-oriented sensor components 
and wherein the orientation of each component is stored with the location of the 
sensor, the polarization analyser comprising: 

means for rotating the observed arrival time signals into a directional 
coordinate system using the stored component orientation information; 

means for user-selection of one of the observed arrival time signals for a trace 
and for providing a timing mark on the trace through which the curve passes through; 

means for retrieving from a parameter file a defined window around the 
selected timing mark; 

means for drawing a line in the map window from the location of the sensor 
for each of the observed arrival time signals in the defined window, in a direction 
defined by the amplitude of the arrival time signal in each of the directions defined 
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by the coordinate system, the length of the line being proportional to the amplitudes 
of the arrival time signals; 

means for moving the trial location along the path defined by the lines, such 
that ambiguity in the event location is removed. 

5 13. Apparatus for determining the location of a seismic event in an earth 

formation zone using seismic traces corresponding to the event collected from one or 
more sensors in the zone, the apparatus being stored on a computer-readable media 
for operation by a user on a computer having a graphical user interface, the apparatus 
comprising: 

10 means for displaying a graphical data window of the seismic traces 

corresponding to the event, the traces representing observed arrival time signals from 
the sensors plotted with respect to time; 

means for displaying a window of a base map of the positions of the one or 
more sensors within the zone; 

15 means for displaying in the graphical data window a curve plotting calculated 

arrival time signals for the traces, the curve being developed from a user-selected one 
of the observed arrival time signals displayed in the graphical data window and from 
a user-selected trial location of the event selected in the base map window; and 

means for user-repositioning of the trial location in the base map window such 

20 that the displayed curve moves in correspondence with the repositioning of the trial 
location, the trial location being repositioned by the user until an acceptable visual 
coincidence of the curve is reached with the observed arrival time signals of the traces 
to thereby yield the location of the event. 

14. The apparatus of claim 13 further including means for refining the 
25 location of the event, comprising: 

means for segmenting the curve, each such segment corresponding to one of 
the traces and representing the calculated arrival time signal for the corresponding 
trace; 

means for user-repositioning of user-selected ones of the curve segments; and 
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means for performing a functional analysis on the calculated arrival time 
signals following user-repositioning to yield a refined location of the event. 

15. The apparatus of claim 14 wherein the means for refining the location 
of the event further comprise: 

means for displaying a set of new curve segments on the traces in the graphical 
display window, the new curve segments representing recalculated arrival time signals 
following performance of the functional analysis; 

means for user-repositioning of user-selected ones of the new curve segments, 
such that a subsequent functional analysis is performed on the calculated arrival time 
signals following the user-repositioning to yield a further refined location of the event. 

16. The apparatus of claim 13, 14 or 15, further including a polarization 
analyser to remove ambiguity in the location of the seismic event, the one or more 
sensors for collecting the event data including differently-oriented sensor components 
and wherein the orientation of each component is stored with the location of the 
sensor, the polarization analyser comprising: 

means for rotating the observed arrival time signals into a directional 
coordinate system using the stored component orientation information; 

means for user-selection of one of the observed arrival time signals for a trace 
and for providing a timing mark on the trace through which the curve passes through; 

means for retrieving from a parameter file a defined window around the 
selected timing mark; 

means for drawing a line in the map window from the location of the sensor 
for each of the observed arrival time signals in the defined window, in a direction 
defined by the amplitude of the arrival time signal in each of the directions defined 
by the coordinate system, the length of the line being proportional to the amplitudes 
of the arrival time signals; 

means for moving the trial location along the path defined by the lines, such 
that ambiguity in the event location is removed. 
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17. A computer apparatus for determining the location of a seismic event 
in an earth formation zone using seismic traces corresponding to the event collected 
from one or more sensors in the zone, the apparatus being stored on a computer- 
readable media for operation by a user on a computer having a graphical user 
5 interface, the apparatus comprising: 

means for displaying a graphical data window of the seismic traces 
corresponding to the event, the traces representing observed arrival time signals from 
the sensors plotted with respect to time; 

means for displaying a window of a base map of the positions of the one or 

10 more sensors within the zone; 

means for generating a map-locating the event by displaying in the graphical 
data window a curve plotting calculated arrival time signals for the traces, the curve 
being developed from a user-selected one of the observed arrival time signals 
displayed in the graphical data window and from a user-selected trial location of the 

15 event selected in the base map window, followed by user-repositioning of the trial 
location in the base map window such that the displayed curve moves in 
correspondence with the repositioning of the trial location, the trial location being 
repositioned by the user until an acceptable visual coincidence of the curve is reached 
with the observed arrival time signals of the traces to thereby yield the location of the 

20 event; and 

means for pick-locating the event to further refine its location after the map- 
locating, including displaying timing marks at positions of the curve on the traces, 
repositioning picked ones of the timing marks to more accurately visually coincide 
with the observed arrival time signals, and performing a functional analysis on the 
25 calculated arrival time signals following repositioning of the timing marks to yield a 
refined location of the event based on the changes in the calculated arrival time 
signals caused by the repositioning of the timing marks. 

18. A method for determining the location of a seismic event in an earth 
formation zone using seismic traces corresponding to the event collected from one or 
30 more sensors in the zone, the method for performance by a user on a computer having 
a graphical user interface, the method comprising: 
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displaying a graphical data window of the seismic traces corresponding to the 
event, the traces representing observed arrival time signals from the sensors plotted 
with respect to time; 

displaying a window of a base map of the positions of the one or more sensors 
within the zone; and 

map-locating the event by displaying in the graphical data window a curve 
plotting calculated arrival time signals for the traces, the curve being developed from 
a user-selected trial location of the event selected in the base map window, the curve 
then being repositioned to visually coincide with the observed arrival time signals of 
the traces and thereby indicate the location of the event based on the changes in the 
calculated arrival time signals caused by the repositioning of the curve. 

19. The method of claim 18 further comprising repositioning the curve 
wherein the curve is graphically moved in correspondence with movement of the trial 
location in the map window. 

20. The method of claim 18 further comprising performing the repositioning 
with an instruction pointer operated by the user. 

21. The method of claim 18 further including pick-locating the event to 
further refine its location after the map-locating, comprising: 

repositioning timing marks for traces of the curve corresponding to the 
calculated arrival time signals at picked ones of the traces, the timing marks being 
repositioned to more accurately visually coincide with the observed arrival time 
signals, and 

performing a functional analysis on the calculated arrival time signals 
following repositioning of the timing marks to yield a refined location of the event 
based on the changes in the calculated arrival time signals caused by the repositioning 
of the timing marks. 

22. The method of claim 21 wherein the functional analysis is a least- 
squares statistical analysis. 
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23. The method of claim 21 further including automatically picking the 
traces for which timing marks are to be repositioned. 

24. The method of claim 18 further including providing options in 
formatting of the display of the seismic traces in the graphical data window. 

5 25. The method of claim 24 wherein the formatting options include one or 

more of plot zooming, changing the trace gain, changing the time scale, and changing 
the spacing between traces. 

26. The method of claim 18 wherein the map locating method further 
comprises: 

10 selecting a trace in the graphical data window wherein the trace becomes 

highlighted; 

selecting an observed arrival time of a trace to provide a timing mark through 
which an analytically predicted time of the curve passes through, such that all 
calculated arrival times of the curve are determined relative to the timing mark; 
15 selecting a trial location of the seismic event in the map window wherein the 

selected trial location is highlighted; 

displaying the coordinates of the location of the event as determined by the 
curve repositioning; and 

storing the calculated arrival times and the event location as determined by the 
20 curve repositioning. 

27. The method of claim 21 wherein the pick-locating method further 
comprises: 

converting the curve to timing marks corresponding to the calculated arrival 
time signals; and 

25 highlighting a trace for which repositioning is to occur. 

28. The method of claim 27 wherein the pick-locating method further 
comprises generating a tabulation of the extent of agreement between the repositioned 
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calculated arrival times and the observed arrival times as indicated by the functional 
analysis performed. 

29. The method of claim 18 further including a polarization analysis to 
remove ambiguity in the location of the seismic event, the one or more sensors for 
collecting the event data including differently-oriented components at the same 
location and wherein the orientation of each component is stored with the location of 
the sensor, the polarization analysis comprising: 

rotating the observed arrival time signals into a directional coordinate system 
using the stored component orientation information; 

user-selection of at least one of the observed arrival time signals for a trace 
and providing a timing mark on the trace through which the curve passes through; 

drawing a line in the map window from the location of the sensor for each of 
the selected observed arrival time signals, in a direction defined by the amplitude of 
the arrival time signal in each of the directions defined by the coordinate system, the 
length of the line being proportional to the amplitudes of the arrival time signals; and 

moving the trial location along the path defined by the lines, such that 
ambiguity in the event location is removed. 

30. A system for determining the location of an extension of a hydraulic 
fracture within an earth formation zone of interest comprising: 

seismic signal sensors placed in an array with respect to the formation zone; 

means for transmitting signals relating to an event caused by extension of the 
fracture from the sensors at a predetermined sampling rate to a memory device; 

means for selectively retrieving the signals from the memory device and 
determining the location of the event by displaying a graphical data window of the 
seismic traces corresponding to the event, the traces representing observed arrival time 
signals from the sensors plotted with respect to time; displaying a window of a base 
map of the positions of the sensors within the zone; and displaying in a graphical data 
window a curve plotting calculated arrival time signals for the traces, the curve being 
developed from a user-selected trial location of the event selected in the base map 
window, the curve then being repositioned to visually coincide with the observed 
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arrival time signals of the traces and thereby indicate the location of the event based 
on the changes in the calculated arrival time signals caused by the repositioning of the 
curve. 

31. A computer program for determining the location of a seismic event in 
5 an earth formation zone using seismic traces corresponding to the event collected from 
one or more sensors in the zone, the program being stored on a computer-readable 
media for operation by a user on a computer having a graphical user interface, the 
program comprising: 

instructions for displaying a graphical data window of the seismic traces 
10 corresponding to the event, the traces representing observed arrival time signals from 
the sensors plotted with respect to time; 

instructions for displaying a window of a base map of the positions of the one 
or more sensors within the zone; and 

instructions for map-locating the event by displaying in the graphical data 
15 window a curve plotting calculated arrival time signals for the traces, the curve being 
developed from a user-selected trial location of the event selected in the base map 
window, the curve then being repositioned to visually coincide with the observed 
arrival time signals of the traces and thereby indicate the location of the event based 
on the changes in the calculated arrival time signals caused by the repositioning of the 
20 curve. 
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